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IV 


ABSTRACT 


A  TOPAZ-II  Thermionic  Fuel  Element  (TFE)  test  stand  is  under  study  at  the 
TOPAZ  International  Program.  Extensive  cross-sectional  drawings  of  this  system  are 
presented  along  with  detailed  descriptions  of  its  construction  A  thermal  analysis  of  the 
test  stand  is  developed  using  a  thermal  resistance  network  method  and  data  collected  on 
the  test  stand.  Critical  resistances  in  the  network  include  the  0.5  mm  regulated  helium 
gap,  cooling  water  channel,  and  cesium  interelectrode  gap.  The  thermal  analysis  indicates 
that  for  an  input  power  to  the  active  zone  of  3003  W,  and  a  regulated  helium  gap  pressure 
of  10  torr,  the  thermal  conductivity  of  the  regulated  gap  helium  ranges  from  0.073-0.087 
W/mK.  This  is  much  less  than  the  thermal  conductivity  calculated  as  a  function  of 
temperature  only  of  0.276  W/mK.  For  a  water  mass  flow  rate  of  0  04  kg/s,  the  convective 
heat  transfer  coefficient  ranges  from  528-1012  W/nrK  based  on  correlations  for  fully 
developed  laminar  and  turbulent  flow.  Our  model  may  be  used  in  conjunction  with  a 
thermionic  emission  model  in  order  to  more  completely  evaluate  the  operating 
characteristics  of  the  test  stand,  and  to  better  correlate  results  with  actual  operating 
conditions  in  the  TOPAZ-II  reactor  core. 
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I.  INTRODUCTION 


The  TOPAZ-II  Thermionic  Fuel  Element  (TFE)  test  stand  was  originally  built  and 
tested  at  the  Scientific  Industrial  Association  "LUTCH"  in  Podolsk,  Russia.  In  December, 
1992,  the  stand  was  reassembled  at  the  New  Mexico  Engineering  Research  Institute 
(NMERI)  in  Albuquerque,  New  Mexico.  The  TFE  test  stand  was  designed  to  conduct 
single  TFE  non-nuclear  testing  by  replacing  the  nuclear  fuel  with  a  tungsten  heating 
element  (TISA  heater).  The  test  stand  provides  a  means  to  test  a  single  TFE  so  that  its 
performance  can  be  better  observed  without  the  interference  and  complexity  of  the 
integrated  TOPAZ-II  system.  This  removes  the  direct  coupling  of  the  TOPAZ-II  in-core 
TFE  with  the  coolant  loop,  radiator  and  the  other  TFEs  electrically  connected  in  series. 
Acceptance  and  demonstration  testing  were  conducted  throughout  1993  (Luchau  et  al., 
1993). 

In  order  to  better  understand  and  predict  the  behavior  of  a  TFE  in  the  test  stand,  a 
detailed  thermal  analysis  based  on  a  network  method  needs  to  be  developed  to  account  for 
the  thermal  losses  in  the  stand.  Because  of  the  large  amounts  of  Russian  technical  data 
and  engineering  drawings  associated  with  the  TOPAZ-II  TFE  test  stand,  and  due  to  the 
extreme  difficulty  in  translating  Russian  technical  information,  the  implementation  of  this 
thermal  model  required  the  creation  of  detailed  and  extensive  cross-sectional  drawings  for 
this  thesis.  Mathcad Plus  5.0  (MathSoft,  Inc.,  1994)  was  the  primary  technical  calculation 
software  used  to  solve  the  thermal  analysis  model. 

Since  the  test  stand  provides  a  means  for  testing  a  single  TFE  without  the  support 
of  the  supporting  TOPAZ-II  reactor  systems,  there  are  several  differences  in  the 
construction  and  operation  of  the  test  stand  from  the  TOPAZ-II  reactor  core.  In  the  test 
stand,  water  cooling  is  used  to  remove  the  excess  heat  instead  of  the  liquid  metal  cooling 
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that  is  used  in  the  reactor  core.  The  cooling  water  flows  through  the  test  stand  water 
jacket  along  a  torturous  path  to  aid  in  its  heat  removal  ability  The  heat  transfer 
coefficient  for  the  cooling  water  is  directly  dependent  on  the  flow  cross-sectional  area  and 
the  total  heat  transfer  area,  both  of  which  are  challenging  to  determine.  Thermocouples 
internal  to  the  test  stand  provide  a  axial  temperature  distribution  along  the  TFE  working 
section.  No  internal  temperature  measurements  are  made  in  the  TOPAZ-II  reactor  core 
TFEs.  Heaters  are  also  located  internally  in  the  test  stand  to  provide  a  means  of  degassing 
the  test  stand  components  prior  to  operation.  A  0.5  mm  gap  in  the  test  stand  just  inside  of 
the  cooling  water  jacket  is  filled  with  low  pressure  helium  to  provide  a  means  of 
controlling  the  heat  transfer  through  the  test  stand.  This  gap  is  used  to  help  the  test  stand 
better  simulate  the  operating  conditions  experienced  by  a  TFE  in  the  reactor  core.  Due  to 
the  low  pressure  band  that  the  gap  is  adjusted  over  (from  1  to  10  torr),  the  pressure 
dependence  of  the  rarefied  helium's  thermal  conductivity  is  difficult  to  define. 

Experimental  data  for  this  work  were  collected  in  March  1994,  with  follow-on 
experimental  data  taken  in  August  1994  at  NMERI.  Along  with  aiding  in  the  development 
of  the  thermal  model  for  the  test  stand,  these  data  are  used  to  investigate  several  of  the 
critical  resistances  from  the  thermal  network.  A  good  understanding  of  the  test  stand 
thermal  characteristics  is  essential  in  evaluating  the  performance  of  a  single  TFE  operating 
in  the  test  stand,  and  correlating  that  performance  to  a  TFE  operating  in  the  TOPAZ-II 
reactor  core.  The  experimentally  measured  input  power  to  the  TFE  working  section  along 
with  the  collector  sleeve  temperature  information  is  also  used  as  an  input  to  the  TOPAZ-II 
Thermionic  Transient  Analysis  Model  (TITAM).  TITAM  is  a  computer  simulation 
program  developed  by  the  University  of  New  Mexico  to  look  at  both  steady-state  and 
transient  operating  characteristics  of  the  TOPAZ-II  reactor  system  (Paramonov  and 
El-Genk,  1994). 


The  goals  of  this  work  include  the  creation  of  detailed  engineering  drawings  for 
the  test  stand,  investigation  of  test  stand  material  properties,  development  of  a  thermal 
model  for  the  test  stand,  and  benchmarking  the  model  with  the  experimental  data 
collected.  Additionally,  experimental  data  are  used  to  empirically  fix  a  range  for  several  of 
the  thermal  properties  that  cannot  be  easily  determined 

LT  J.  R.  Venable  is  writing  a  companion  thesis  that  will  concentrate  on  the 
thermionic  theory  behind  the  data  collected  for  this  thesis.  Venable's  thesis  will  be 
available  after  March  of  1995  and  a  joint  paper  will  be  presented  ("Operational  Testing 
and  Thermal  Modeling  of  a  TOPAZ-II  Single-Cell  Thermionic  Fuel  Element  Test  Stand") 
at  the  Space  Nuclear  Power  and  Propulsion  Symposium  held  in  Albuquerque,  New 
Mexico,  on  January  8-12,  1995. 


II.  BACKGROUND 


The  continued  exploration  and  exploitation  of  space  will  require  increasing 
amounts  of  power  as  technology  continues  to  advance  Supporting  high  definition 
television  broadcasts  and  the  rapidly  growing  world-wide  communication  needs,  moving 
large  level  research  and  industrial  opportunities  to  Earth  orbit  or  into  the  cislunar  region, 
revisiting  and  colonizing  the  Moon,  and  pushing  the  envelope  of  civilization  to  Mars  and 
beyond  are  future  goals  which  will  demand  large  electrical  power  sources  that  are  both 
compact  and  have  long  lifetimes  These  larger  power  sources  will  need  to  survive  harsh 
environments,  including  radiation  belts,  the  vacuum  of  deep  space,  and  hostile 
atmospheres  of  other  planets  and  natural  satellites  in  the  solar  system.  Nuclear  power  can 
provide  the  compact  power  sources  necessary  to  meet  these  future  needs. 

A.         UNITED  STATES  SPACE  NUCLEAR  POWER  PROGRAMS 

Since  the  mid-1950s,  the  United  States  has  been  interested  in  the  use  of  nuclear 
power  in  space  applications.  From  radioisotope  power  systems,  nuclear  rockets,  and  the 
use  of  nuclear  reactors  for  power  production,  interest  in  the  space  nuclear  power  area  was 
motivated  by  a  favorable  political  environment  and  public  support  until  the  early  1970s. 

The  Transit  4 A  navigational  satellite  was  launched  in  June  of  1961  with  a 
SNAP-3B  Radioisotope  Thermoelectric  Generator  (RTG)  to  provide  power.  Since  that 
first  successful  launch  of  a  nuclear  power  source,  a  total  of  24  launches  using  RTGs  have 
followed  including  the  Pioneer,  Viking,  and  Voyager  missions.  Although  the  RTG  has 
proved  to  be  very  useful  for  many  applications,  the  ability  to  support  a  very  ambitious 
space  program  required  the  development  of  a  nuclear  reactor  based  on  the  principle  of 
atomic  fission  for  both  propulsion  purposes  and  to  provided  electrical  power. 


Project  Rover  was  initiated  at  Los  Alamos  National  Laboratory  in  1955  and  was 
the  beginning  of  approximately  18  years  of  research  and  development  into  nuclear  rockets. 
Although  the  nuclear  rocket  program  was  terminated  in  late  1972,  it  was  the  start  of  the 
search  for  a  high  temperature  nuclear  reactor.  The  research  and  demonstration  test  of  the 
Nuclear  Engine  for  Rocket  Vehicle  Application  (NERVA)  program  proved  that  a  nuclear 
reactor  could  provide  a  high  power  density  rocket  with  high  thrust  and  specific  impulse. 
(Angelo  and  Buden,  1985) 

Active  development  of  a  nuclear  reactor  to  provide  electric  power  began  in  1957 
with  the  Systems  for  Nuclear  Auxiliary  Power  (SNAP)  program.  These  nuclear  reactors 
required  a  power  generation  method  to  convert  the  heat  created  by  atomic  fission  into 
electrical  power.  The  initial  concept  was  SNAP-2,  a  thermal  uranium  zirconium-hydride 
reactor  that  utilized  a  Rankine  cycle  with  mercury  as  the  working  fluid  to  convert  the 
nuclear  heat  to  electricity.  This  was  the  precursor  to  two  separate  concepts,  the  SNAP-8 
and  the  SNAP- 10 A,  both  of  which  were  initiated  in  1960.  The  SNAP-8  reactor  utilized  a 
mercury  Rankine  cycle  along  with  a  larger  version  of  the  SNAP-2  reactor  system  to 
support  anticipated  higher  power  demands.  (Angelo  and  Buden,  1985) 

The  SNAP-10A  was  also  based  on  the  SNAP-2  reactor  system,  but  instead  of 
using  the  dynamic  Rankine  power  conversion  cycle,  it  used  a  thermoelectric  power 
conversion  system.  The  use  of  static  power  conversion  provides  a  means  of  converting 
heat  directly  into  electricity  without  the  need  of  a  secondary  working  fluid  (such  as 
mercury  for  a  Rankine  cycle  or  helium-xenon  for  a  closed  Brayton  cycle)  and  the 
associated  turbines  and  pumps/compressors.  In  April  of  1965,  the  SNAPSHOT  mission 
was  successfully  launched  carrying  a  SNAP-10A  reactor  to  approximately  a  1300  km 
orbit.  Although  the  SNAP-10A  reactor  shutdown  only  43  days  after  becoming  critical,  it 
proved  that  a  nuclear  reactor  could  be  successfully  launched  and  operated  on  orbit.   This 


was  the  only  attempt  by  the  United  States  to  launch  a  nuclear  reactor  into  space    (Angelo 
and  Buden,  1985) 

Another  static  power  conversion  method  that  was  developed  in  parallel  with  the 
thermoelectric  method  was  thermionic  power  conversion  Whereas  thermoelectrics  work 
on  the  basis  of  a  thermally  induced  emf  at  the  junction  of  two  dissimilar  semiconductors, 
thermionics  works  basically  by  "boiling  off'  electrons  from  a  hot  emitter  surface  and 
collecting  the  electrons  on  a  relatively  cooler  collector  surface  across  a  small  gap.  The 
collector  and  emitter  are  then  electrically  connected  together  through  an  external  load 
which  completes  the  electrical  circuit.  The  process  of  thermionics  is  discussed  in  more 
detail  in  the  next  chapter.  Although  thermionic  devices  were  tested  in  the  mid-  to 
late- 1960s,  no  complete  reactor  systems  using  thermionics  was  actually  built  in  the  United 
States.  (Angelo  and  Buden,  1985) 

Due  to  a  changing  political  environment  and  shifting  emphasis  in  the  space 
program,  all  nuclear  reactor  development  programs  were  canceled  at  the  end  of  1972. 
Although  some  research  did  continue  after  1972,  it  wasn't  until  the  early  1980s  that 
increased  interest  was  returned  to  space  nuclear  reactors.  This  renewed  interest  was 
mainly  due  to  the  need  for  a  large  power  source  to  support  the  Strategic  Defense  Initiative 
(SDI).  The  advance  space  nuclear  power  program  was  reinitiated  as  the  Space  Power 
Advanced  Reactor  (SPAR)  system  later  renamed  SP-100.  Several  reactor  designs  and 
power  conversion  systems  were  initially  considered  for  SP-100.  The  primary  concept  that 
SP-100  developed  utilized  a  thermoelectric  power  conversion  system.  (Angelo  and 
Buden,  1985)  Research  on  SP-100  has  continued  until  1994  when  once  again  the  political 
environment  has  changed  and  funding  for  the  program  is  being  terminated. 


B.         FORMER  SOVIET  UNION  SPACE  NUCLEAR  POWER  PROGRAMS 

The  former  Soviet  Union,  like  the  United  States,  was  very  active  in  the  research 
and  development  of  space  nuclear  power.  From  the  early  1960s  to  1989,  a  total  of  38 
nuclear  reactor  systems  were  launched  by  the  former  Soviet  Union.  The  majority  of  the 
nuclear  reactors  launched  were  for  the  Radar  Ocean  Reconnaissance  SATellites 
(RORSATs)  and  utilized  thermoelectric  energy  conversion.  The  other  two  nuclear 
reactors  launched  were  COSMOS-1818  and  COSMOS- 1867.  Both  COSMOS  satellites 
were  test  platforms  for  the  TOPAZ-I  reactor  system  utilizing  multi-cell  thermionic  energy 
conversion.  (Polansky  et  al.,  1994) 

In  the  early  1960s,  the  former  Soviet  Union  began  development  of  single-cell  and 
multi-cell  TFEs.  In  1967  the  two  concepts  were  sent  for  system  development  to  two 
separate  technical  institutes.  The  single-cell  TEE  was  developed  by  the  Central  Design 
Bureau  of  Machine  Building  (CDBMB)  in  St.  Petersburg  under  the  project  name  of 
ENISY.  The  multi-cell  TFE  was  developed  by  the  Scientific  and  Production  Cooperation 
"Krasnaya  Zvezda"  in  Moscow  under  the  project  name  of  TOPAZ,  the  Russian  acronym 
for  Thermionic  Experiment  with  Conversion  in  Active  Zone.  Both  concepts  were 
developed  cooperatively  between  the  two  institutes  up  until  the  early  1 970s.  This  lead  to 
many  similarities  between  the  two  reactor  designs.  Cooperative  efforts  ended  in  the  early 
1970s  due  to  differences  between  the  two  institutes  and  the  programs  developed 
independently.  (Voss,  1994) 

In  1987  and  1988,  the  multi-cell  TFE  reactor  TOPAZ  was  successfully  launched 
twice  on  COSMOS-1818  and  COSMOS- 1867.  The  first  mission  operated  for  142  days 
and  the  second  mission  operated  for  342  days  before  the  planned  depletion  of  the  cesium 
supplies  on  board  terminated  the  tests.  Both  test  were  extremely  successful  and  confirmed 
that  the  on-orbit  data  correlated  well  to  ground  test  results.    (Paramonov  and  El-Genk, 


1994)  Although  there  were  26  single-cell  TFE  ENISY  reactor  systems  manufactured  and 
19  systems  tested  between  1970  and  1989,  none  of  the  single-cell  TFE  technology  were 
ever  launched  and  tested  on  orbit.  (Schmidt  et  al.,  1994) 

C.         TOPAZ  INTERNATIONAL  PROGRAM 

At  the  end  of  the  1980s,  the  Russians  approached  a  private  U.S.  company 
(International  Scientific  Products,  San  Jose,  CA)  about  transferring  Russian  thermionic 
technology  through  the  purchase  of  Russian  space  reactor  hardware.  During  the 
negotiations,  the  U.S.  company  inadvertently  renamed  ENISY  to  TOPAZ-II  and  TOPAZ 
to  TOPAZ-I  (Voss,  1994).  In  May  of  1992,  the  U.S.  Government  purchased  two 
complete  TOPAZ-II  reactor  systems  and  associated  test  and  support  hardware  through 
the  purchasing  agent  International  Scientific  Products.  Funding  was  provided  by  the 
Strategic  Defense  Initiative  Organization  (SDIO)  to  investigate  the  possible  use  of  the 
TOPAZ-II  reactor  system  in  support  of  the  Nuclear  Electric  Propulsion  (NEP)  Space  Test 
Mission.  This  project  has  since  been  named  the  Topaz  International  Program  (TIP). 
(Wyant,  1994) 

The  TIP  facility  is  located  at  NMERI  which  is  owned  and  operated  by  the 
University  of  New  Mexico  in  Albuquerque,  NM.  Funding  is  still  being  provided  by  the 
former  SDIO  (now  the  Ballistic  Missile  Defense  Organization,  BMDO)  and  administered 
by  the  U.S.  Air  Force's  Phillips  Laboratory.  The  New  Mexico  alliance  for  Thermionics 
which  includes  the  University  of  New  Mexico,  Sandia  National  Laboratories  and  Los 
Alamos  National  Laboratory  provides  research  support  for  the  TIP  program.  The 
international  effort  includes  scientist  and  engineers  from  the  Russia,  England,  and  France. 
Although  the  flight  program  for  TIP  has  been  canceled,  in  March  of  1994  four  additional 
TOPAZ-II  units  were  purchased  for  continuing  research  into  thermionics  and  reactor 
operation.  (Wyant,  1994) 
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III.  THERMIONIC  ENERGY  CONVERSION 


A. 


BASIC  THEORY 


The  thermionic  process  can  be  considered  to  be  a  "heat  engine"  which  utilizes 
electrons  as  its  working  fluid  to  convert  heat  into  electricity.  The  thermionic  converter 
consists  of  a  emitter,  collector,  and  a  small  gap  between  the  two.  The  collector  is 
electrically  connected  to  the  emitter  through  an  external  load  to  complete  the  electrical 
circuit.  Figure  3-1  shows  a  simplified  flat-plate  thermionic  converter  (Rasor,  1991). 
Electrons  are  "boiled"  off  the  hot  emitter  and  travel  across  the  gap  and  are  "condensed"  on 
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Figure  3-1 .  Simplified  Flat-plate  Thermionic. 
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the  relatively  cooler  collector.  The  potential  difference  between  the  emitter  and  collector 
creates  a  driving  force  which  causes  a  current  to  flow  through  the  external  load.  (Angelo 
andBuden,  1985) 

Since  the  thermionic  device  acts  as  a  "heat  engine,"  the  maximum  efficiency  of  the 
thermionic  cycle  can  be  characterized  as  the  cycles  Carnot  efficiency.  Some  of  the 
inefficiencies  of  the  thermionic  cycle  include:  (1)  radiation  heat  transfer  from  the  emitter 
to  the  collector;  (2)  negative  space  charge  in  the  gap  that  interferes  with  the  electron  flow 
from  the  emitter  to  the  collector;  and  (3)  thermal  losses.  These  inefficiencies  prevent  the 
thermionic  cycle  from  approaching  Carnot  cycle  efficiency.  (Angelo  and  Buden,  1985) 

Radiation  heat  transfer  from  the  emitter  to  the  collector  cannot  be  mitigated. 
Higher  collector  temperatures  would  not  only  reduced  the  amount  of  radiation  heat 
transfer  but  would  also  allow  for  a  much  smaller  radiator  since  the  radiators  surface  area  is 
inversely  proportional  to  the  radiators  operating  temperature  to  the  fourth  power. 
However,  as  the  collector  temperature  increases,  the  number  of  electrons  being  "boiled" 
off  the  collectors  surface,  called  electron  back  emission,  also  increases.  This  in  effect 
causes  the  overall  cycle  efficiency  to  decrease  dramatically.  A  design  trade-off  must  be 
made  to  understand  the  optimum  collector  operating  temperature.  (Angelo  and  Buden, 
1985) 

The  negative  space  charge  in  the  gap  can  be  decreased  by  either  making  the  gap 

size  smaller  or  by  injecting  a  low  pressure  gas  into  the  gap.     A  vacuum  diode  is  a 

thermionic  device  that  relies  on  a  very  small  gap  between  the  electrodes  to  reduce  the 

current  inhibiting  effects  of  the  electric  field  produced  by  the  negative  space  charge. 

Theoretical  spacing  needed  to  produce  useful  current  density  approaches  the  associated 
Debye  length  (~  1  urn)  which  is  for  a  practical  purpose  very  difficult  to  manufacture  and 
maintain.    This  very  narrow  gap  can  be  widened  somewhat  (up  to  at  most  0.5  mm)  by 

injecting  a  gas  plasma  into  the  gap.   This  type  of  thermionic  device  is  known  as  a  plasma 
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diode.  The  primary  purpose  of  this  plasma  is  to  facilitate  to  transfer  of  electrons  through 
the  gap  while  minimizing  the  potential  drop  and  current  attenuation  due  to  the  presence  of 
the  plasma.  (Rasor,  1991) 

The  most  common  gas  that  is  used  in  this  type  of  diode  is  cesium  which  has  the 
lowest  ionization  potential  of  any  chemical  element  (Angelo  and  Buden,  1985).  A 
measure  of  the  amount  of  energy  necessary  for  a  electron  to  escape  through  the  potential 
barrier  associated  with  the  surface  of  a  metal  is  called  the  work  function  and  is  given  in 
units  of  electron-volts  (eV)  (Angrist,  1976).  Not  only  does  the  cesium  plasma  in  the  gap 
decrease  the  space  charge,  it  also  may  have  the  effect  of  forming  a  monolayer  of  cesium 
on  both  the  emitter  and  collector  causing  a  decrease  in  their  associated  work  functions. 
This  decrease  in  the  metals  work  function  is  directly  related  to  the  pressure  of  the  cesium 
in  the  gap  and  the  electrode  temperature.  Cesium  pressure  can  be  controlled  by  adjusting 
the  temperature  of  the  liquid  cesium  reservoir  which  in  effect  controls  the  vapor  pressure 
of  the  evaporating  cesium.  (Angelo  and  Buden,  1985) 

B.         TOPAZ-II  NUCLEAR  REACTOR 

The  TOPAZ-II  nuclear  reactor  consists  of  37  individual  TFEs  incorporated  in  a 
cylindrical  core  arrangement.  Figure  3-2  shows  a  cross-sectional  view  of  the  TOPAZ-II 
reactor.  Each  of  37  TFEs  consist  of  a  UO,  fuel  assemble,  a  single-cell  cylindrical 
thermionic  device,  and  a  sodium-potassium  (NaK)  liquid  metal  cooling  channel.  Electrical 
power  is  supplied  to  the  loads  by  34  of  the  TFEs  connected  in  series.  These  TFEs  are 
designed  to  produce  5-6  kWe  at  27  VDC  for  up  to  3  years.  The  other  three  TFEs  are 
connected  in  parallel  to  provide  power  to  the  electromagnetic  pump.  This  pump  is  used  to 
generate  the  pumping  power  necessary  to  circulate  the  NaK.  The  overall  length  of  the 
reactor  core  is  37.5  cm  with  an  overall  width  of  26.0  cm.  (Paramonov  and  El-Genk, 
1994) 
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Figure  3-2.  Cross-sectional  View  of  the  TOPAZ-II. 

The  TFEs  are  inserted  into  axial  channels  in  a  zirconium-hydride  (ZrH)  moderator. 
Twelve  reflector  drums  surround  the  core,  three  are  safety  drums  which  are  turned  out 
prior  to  initial  criticality  when  on  orbit.  The  other  nine  reflector  drums  are  used  to  control 
the  nuclear  reaction  during  operation  and  are  called  control  drums.  All  of  the  reflector 
drums  are  made  out  of  beryllium  clad  in  stainless  steel  with  a  small  section  of  neutron 
absorbing  material  made  of  boron-carbide  (B4C).  A  beryllium  reflector  surrounds  the 
remainder  of  the  core  as  shown  in  Figure  3-2.  (Schmidt  et  al.,  1994) 


Excess  heat  is  removed  from  the  TFEs  using  a  eutectic  NaK  coolant  pumped 
through  stainless  steel  piping  by  the  electromagnetic  pump.  The  hot  coolant  flows 
through  a  conical  radiator  where  the  waste  heat  is  rejected  to  space.  During  normal 
operation  the  thermionic  emitter  temperature  ranges  between  1800  K  to  2100  K  and  the 
NaK  coolant  maintains  the  collector  temperature  between  743  K  to  843  K.  (Schmidt  et 
al.,  1994) 
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c. 


TOPAZ-II  THERMIONIC  FUEL  ELEMENT 


The  TOPAZ-II  TFE  consists  of  a  single-cell  cylindrical  thermionic  converter.  A 
single-cell  converter  utilizes  the  entire  length  of  the  reactor  core  as  a  single  device.  This  is 
in  contrast  to  the  TOPAZ-I  TFE  which  is  made  up  of  multi-cell  thermionic  converters.  A 
simple  comparison  of  the  electrical  circuit  for  both  a  single-cell  and  multi-cell  TFE  are 
shown  in  Figure  3-3.  The  advantages  of  using  multi-cell  converters  are  that  they  produce  a 
higher  power  density  and  higher  conversion  efficiency.  This  is  mainly  due  to  less  ohmic 
losses  in  the  multi-cell  because  of  the  shorter  thermionic  elements.  Another  advantage  of 
the  multi-cell  design  is  that  each  of  the  cells  are  connected  in  series  resulting  in  a  TFE 
output  voltage  in  the  volt  range.  This  reduces  the  amount  of  electrical  losses  in  the  leads 
and  makes  the  power  distribution  system  more  efficient.  (Paramonov  and  El-Genk,  1994) 
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Figure  3-3.  Simplified  Electrical  Schematics  for  Multi-Cell  and  Single-Cell  TFEs. 
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The  two  major  disadvantages  of  using  a  multi-cell  TFE  are  that  manufacturing  is 
significantly  more  complicated  and  the  design  does  not  provide  a  means  of  removing  the 
nuclear  fuel  and  performing  non-nuclear  testing  Although  the  single-cell  TFE  design  has 
a  much  lower  efficiency,  it  does  provide  a  less  complicated  converter  and  a  means  of 
performing  non-nuclear  testing.  This  testing  is  accomplished  by  replacing  the  UO,  fuel 
assembly  with  a  tungsten  heater  (TISA  heater).  (Paramonov  and  El-Genk,  1994)  The 
TISA  heaters  provide  a  temperature  distribution  similar  to  that  created  using  nuclear 
heating.  Detailed  information  on  the  use  of  electrical  heaters  to  simulate  nuclear  heating 
can  be  found  in  El-Genk  and  Xue,  1994. 

The  interelectrode  gap  between  the  emitter  and  collector  is  filled  with  cesium 
vapor  during  normal  operation.  The  cesium  is  supplied  by  a  single  liquid  cesium  reservoir 
that  is  heated  by  one  of  the  NaK  coolant  return  lines.  The  cesium  in  the  gap  is  constantly 
vented  to  space  at  a  rate  of  approximately  0.5  g/day.  This  is  one  of  the  limiting  factors  in 
the  reactor  lifetime.  (Schmidt  et  al.,  1994)  During  normal  operation,  the  cesium  pressure 
in  the  gap  is  maintained  around  1.5-2.0  torr  (Paramonov  and  El-Genk,  1994). 

Axial  and  radial  temperature  distributions  in  the  reactor  core  cause  performance 
differences  in  the  37  TFEs.  The  exact  performance  characteristic  of  each  individual  TFE 
is  difficult  to  determine  due  to  the  following:  (1)  the  thermal  coupling  of  each  TFE  with 
the  TFEs  surrounding  it,  or  lack  of  TFEs;  (2)  the  electrical  coupling  of  each  TFE  with  the 
other  TFEs,  either  in  series  for  the  34  load  power  producing  TFEs  or  in  parallel  for  the  3 
pump  power  producing  TFEs;  (3)  all  TFEs  share  a  common  heat  removal  system;  and  (4) 
all  TFEs  share  a  common  cesium  system.  To  understand  the  individual  thermionic 
performance  of  a  TFE,  the  TOPAZ-II  TFE  test  stand  must  be  used.  This  test  stand  allows 
detailed  analysis  of  a  single  TFE  under  unique  operating  characteristics. 
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IV.  TEST  STAND  DESCRIPTION 


The  TOPAZ-II  TFE  test  stand  is  used  to  perform  non-nuclear  testing  of  a  single 
TOPAZ-II  TFE.  This  allows  the  TFE  performance  to  be  determined  without  the 
interference  of  the  rest  of  the  TOPAZ-II  reactor  systems  and  without  any  interfering 
effects  caused  by  the  TFEs  relative  location  in  the  core  (see  Figure  3-2).  Interference  by 
the  reactor  systems  include  the  direct  thermal  coupling  of  the  TFE  to  the  NaK  coolant, 
coolant  pump,  and  radiator,  and  the  electrical  coupling  of  the  TFE  to  the  other  TFEs 
connected  in  series  in  the  reactor  core.  The  test  stand  also  allows  varying  certain 
parameters  which  could  not  be  varied  on  the  entire  reactor  system  to  investigate  their 
effect  on  TFE  performance.  These  parameters  include  varying  cesium  vapor  pressure  in 
the  interelectrode  gap,  the  removal  of  all  cesium  vapor  from  the  interelectrode  gap, 
changing  the  collector  sleeve  temperature  (directly  effecting  the  collectors  temperature), 
and  operating  the  TFE  at  its  operational  limits  to  determine  break-down  characteristics. 

The  test  stand  contains  the  test  apparatus  within  a  vacuum  chamber  along  with  the 
associated  vacuum,  cesium,  helium,  and  cooling  water  systems  external  to  the  vacuum 
chamber.  As  with  almost  all  Russian  structural  components,  the  vacuum  chamber  is 
constructed  mainly  of  Russian  stainless  steel  type  12X18H10T.  Phillips  Laboratory  has 
determined  that  this  type  of  Russian  stainless  steel  is  closely  approximated  by  the 
American  Iron  and  Steel  Institute  (AISI)  standard  321.  The  vacuum  system  is  used  to 
initially  evacuate  the  vacuum  chamber,  the  helium  cavities,  the  cesium  vapor  system,  and 
the  TISA  heater  channel.  The  test  stand  start-up  procedure  requires  an  involved 
evacuation  and  out-gassing  procedure.  Out-gassing  is  performed  by  numerous  heaters 
throughout  the  test  stand. 
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The  cesium  system  is  used  to  supply  cesium  vapor  to  the  interelectrode  gap  and  to 
maintain  the  cesium  vapor  pressure.  The  system  also  provides  a  means  to  purify  the 
cesium  inventory.  When  the  test  stand  is  operating,  cesium  vapor  is  provided  to  the 
interelectrode  gap  by  a  cesium  reservoir  that  is  inside  the  vacuum  chamber.  This  isolates 
the  internal  cesium  system  from  the  rest  of  the  system  outside  of  the  vacuum  chamber. 
Cesium  vapor  pressure  in  the  interelectrode  gap  is  regulated  and  varied  by  adjusting  the 
inner  cesium  reservoir  saturation  temperature.  An  external  cesium  reservoir  is  provided  to 
evacuate  and  fill  the  cesium  system  when  the  test  stand  is  not  in  operation.  (Luchau  et  al., 
1993) 

The  helium  system  provides  helium  to  the  cavities  and  gaps  inside  the  test  stand,  a 
means  to  purify  the  helium  supply  and  a  means  to  evacuated  the  helium  from  the  system. 
Helium  is  used  to  fill  in  the  gaps  in  the  test  apparatus  to  facilitate  and  control  the  heat 
transfer  process.  Collector  sleeve  temperatures  can  be  controlled  by  varying  the  pressure 
in  the  regulate  helium  gap.  The  mechanism  by  which  the  thermal  conductivity  of  helium  is 
affected  by  pressure  is  discussed  in  detailed  in  the  next  chapter. 

The  water  cooling  system  picks  up  the  excess  heat  throughout  the  test  stand.  The 
components  cooled  by  the  water  system  include  electric  drives  for  various  valves,  vacuum 
pumps,  cesium  reservoirs,  transistor  load  bank  and  the  test  apparatus.  The  volumetric 
flowrate  of  the  cooling  water  flow  to  the  test  apparatus  is  measured  using  a  turbine  type 
flowmeter  in  the  inlet  line.  The  inlet  and  outlet  cooling  water  temperature  is  measured  by 
type-K,  surface  contact  thermocouples  on  the  inlet  and  outlet  stainless  steel  piping.  Two 
thermocouples  are  used  on  the  outlet  piping,  one  on  the  outlet  pipe  that  is  in  thermal 
contact  with  the  vacuum  chamber  and  the  other  on  the  outlet  pipe  that  is  not  in  thermal 
contact  with  the  vacuum  chamber.  The  experimental  data  show  that  both  of  these 
temperature  measurements  are  well  within  the  thermocouple  tolerance  of  each  other.  The 
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flowrate  and  temperature  measurement  accuracies  are  critical  to  the  thermal  modeling  of 
the  test  stand. 

The  stand  utilizes  an  automatic  control  system,  an  automated  data  acquisition 
system,  and  a  variable  electrical  load  bank.  The  entire  control  and  data  acquisition  system 
was  completely  upgraded  from  the  original  Russian  system  to  a  United  States  technology 
based  system.  This  upgrade  occurred  between  the  first  set  of  tests  in  March  of  1994  and 
the  second  set  in  August  of  1994;  however,  the  testing  parameters  were  maintained  the 
same  during  both  data  sets.  The  load  bank  is  composed  of  both  a  transistive  and  a 
resistive  load  bank.  The  load  bank  can  be  controlled  by  a  micro-computer  to  provide 
current-voltage  (I-V)  sweeps  for  the  TEE.  These  sweeps  are  important  in  determining  the 
actual  operating  characteristics  of  the  TFE.  (Luchau  et  al.,  1993) 

A  diagram  of  the  test  apparatus  is  shown  in  Figure  4-1  (based  on  Luchau  et  al., 
1993).  The  upper  and  lower  helium  chambers  provide  helium  to  the  unregulated  gaps  in 
the  test  section  and  to  the  spaces  in  the  TISA  heater.  The  cesium  for  the  inter-electrode 
gap  enters  the  test  section  through  the  lower  helium  chamber  which  is  heated  during 
testing  to  prevent  "cold-spots"  in  the  cesium  line.  Since  the  lower  helium  chamber  is 
heated  and  the  upper  helium  chamber  is  not  heated,  heat  losses  out  of  the  top  and  bottom 
of  the  test  apparatus  are  different.  This  difference  in  the  end  losses  will  affect  the  test 
stand's  thermal  model.  The  upper  and  lower  helium  chambers  also  contain  the  electrical 
connections  to  both  the  emitter  and  collector.  Electrical  connections  to  the  TISA  heater 
are  made  in  a  chamber  above  the  upper  helium  chamber. 

A  lengthwise  cross-sectional  view  of  the  test  stand's  working  section  is  shown  in 
Figure  4-2.  This  figure  shows  the  relative  proportions  of  the  various  test  section 
concentric  cylinders.  The  overall  length  of  the  thermionic  working  section  is  375.0  mm. 
The  length  of  the  TISA  heater  is  300.0  mm  and  the  length  of  the  cooling  water  jacket  is 


19 


o 

o 

m 
i 

in 
CD 
O 

o' 


ct> 
o 

H 
O 


W 

H 

re 

CO 

I-* 


20 


^  a; 

or-  o  n  d  o  rxi  c 

if)    -i    O    01    CO  f-    CD   ^ 

u 


OJ 

CJ 

u 

C 

0) 

(0 

*J 

-n 

(0 

*j 

* 

OJ 

3 

-^ 

O 

u 

0 

V 

(Q 

it 

3 

—i 

0 

0 

[x. 

', 

< 

OJ 

111 

(J) 

u 

(0 

E- 

u 
U 

n> 

, , 

I 

Id 

u 
0) 

< 

C 
0 

ID 
OJ 

f- 

u 

X 

2      cd 
O 

1) 

cy> 

in 
in 

O 


U 


i 
-f 

CU 


U-. 


21 


approximately  326.0  mm.  The  extremely  small  dimensions  and  close  tolerances  gives  an 
indication  of  the  difficulty  associated  with  the  manufacturing  of  the  test  section.  A  radial 
cut-away  view  of  the  complete  test  stand  is  shown  in  Figure  4-3.  An  enlarged  view  of  the 
portion  of  the  test  stand  from  the  collector  sleeve  inward  is  shown  in  Figure  4-4.  The 
overall  diameter  of  the  entire  test  section  is  only  87.0  mm  giving  an  L/D  of  4.31  for  the 
working  section. 

The  TFE  and  its  associated  TISA  heater  are  fit  into  the  collector  sleeve  tube  of  the 
test  stand  as  shown  in  Figure  4-4.  The  TISA  heater  is  constructed  of  a  special  tungsten 
heating  element  and  flexible  lobes  that  act  as  thermocompensators.  The  maximum  supply 
voltage  is  29  VAC  with  a  maximum  current  of  170  A.  The  TISA  heater  can  supply  a 
maximum  heating  power  to  the  working  section  of  4500  W  (Wold  et  al.,  1993).  The  inner 
conductor  of  the  TISA  heater  has  a  diameter  of  6.5  mm  and  the  outer  conductor  has  a 
thickness  of  0.4  mm  and  an  outer  diameter  of  7.0  mm.  The  length  of  the  heater  is  300.0 
mm  which  makes  the  heater  37.5  mm  shorter  at  each  end  than  the  actual  thermionic 
working  section.  This  shorter  length  simulates  the  end  effects  associated  with  nuclear 
heating.  Detailed  discussions  of  the  correspondance  of  TISA  heating  to  nuclear  fuel 
heating  are  presented  in  El-Genk  and  Xue  (1993). 

The  thermionic  device  consists  of  concentric  emitter  and  collector  cylinders 
separated  by  an  interelectrode  gap  and  is  375.0  mm  in  length.  The  emitter  cylinder  is  1.15 
mm  thick  with  an  outer  diameter  of  19.6  mm.  The  emitter  is  made  of  a  monocrystalline 
molybdenum  -  3%  niobium  alloy  with  its  exterior  surface  coated  with  a  0.1  mm  thick 
monocrystalline  tungsten  layer.  The  0.5  mm  interelectrode  gap  is  filled  with  cesium  vapor 
during  normal  operation.  The  interelectrode  gap  is  maintained  by  the  use  of 
scandium-oxide  (Sc203)  spacers  which  prevent  shorting  between  the  emitter  and  collector. 
A  total  of  six  spacer  sets  are  located  axially  along  the  TFE  working  section.  These 
spacers  are  not  shown  in  any  of  the  figures.  The  collector  cylinder  is  1.4  mm  thick  with  an 
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inner  diameter  of  20  6  mm.  The  collector  is  made  of  polycrystalline  molybdenum  with  its 
outer  surface  coated  by  a  0.15  mm  layer  of  alumina  (AL03)  The  alumina  electrically 
insulates  the  TFE  from  the  remainder  of  the  test  rig.  The  thermionic  device  and  TISA 
heater  used  in  the  test  stand  are  identical  in  construction  to  that  used  in  the  actual 
TOPAZ-II  reactor 

A  0.05  mm  gap  separates  the  alumina  insulator  from  the  collector  sleeve  tube. 
This  gap  allows  the  TFE  to  be  fit  into  the  test  rig  and  is  filled  with  helium  during  normal 
operation.  The  helium  is  provided  by  the  upper  and  lower  helium  chambers  (see  Figure 
4-1).  The  helium  pressure  in  this  gap  is  unregulated  and  is  normally  in  the  200  to  300  ton- 
range. 

The  collector  sleeve  tube  is  3.0  mm  thick  with  an  outer  diameter  of  29.9  mm.  It  is 
constructed  of  Russian  1X18H10T  stainless  steel.  The  collector  sleeve  contains  12  evenly 
spaced,  2.0  mm  deep  slots  in  its  outer  surface  that  run  the  length  of  the  collector  sleeve. 
Figure  4-5  shows  the  relative  location  of  the  thermocouple  slots  in  the  collector  sleeve. 
These  slots  contain  the  only  thermocouples  located  in  the  interior  of  the  test  rig.  The 
thermocouples  are  positioned  along  the  collector  sleeve  to  provide  an  axial  temperature 
distribution  in  the  working  section.  The  thermocouple  wires  are  packed  in  alumina  inside 
of  stainless  steel  tubes  which  are  laid  inside  the  slots.  These  thermocouples  are  positioned 
in  the  test  rig  at  the  location  where  the  NaK  liquid  metal  cooling  flows  through  the  TFE  in 
the  actual  TOPAZ-II  reactor.  This  allows  correlation  of  the  test  stand  data  to  data 
obtained  from  the  complete  reactor  testing.  Another  0.05  mm  helium  gap  separates  the 
collector  sleeve  tube  from  the  copper  conductor.  This  small  gap  is  for  assembly  purposes 
only.  The  helium  in  this  gap  is  also  supplied  from  the  upper  and  lower  helium  chambers 
and  is  unregulated  in  pressure. 


25 


.  v  '  '  v  \  \  v  \  \  \  \  ■  ".'.'^'au  \ \  '  y ■■  v  \  *.  v  v  v  \  v  \  w  \  \  ■.  ■>  ■  ■  v  s  5  5 


^1 


i 


O 
O 

o 

c 

O* 

r 

o 
o 

a 
§ 

S- 

o 

o 

e= 

I 

a* 

1 


26 


The  copper  conductor  is  15  0  mm  thick  with  an  outside  diameter  of  60.0  mm.  It  is 
made  of  99.99%  pure  copper  and  provides  a  heat  transfer  medium  outward  to  the  cooling 
water  The  outer  surface  contains  five  radially  spaced  slots  that  each  contain  four 
nichrome  (NiCr)  heating  elements.  These  heaters  are  used  for  outgassing  the  rig  prior  to 
startup  and  are  secured  during  normal  operation  The  heater  wires  are  insulated  with 
alumina  and  sheathed  in  stainless  steel.  A  0.5  mm  gap  separates  the  copper  conductor 
from  the  inner  cooling  water  jacket.  This  gap  is  filled  with  helium  supplied  from  a  source 
external  to  the  test  rig.  The  pressure  of  the  helium  is  regulated  in  the  band  of  1  to  10  torr 
during  normal  operation  to  control  the  collector  sleeve  temperature  by  varying  the  thermal 
resistance  of  the  gap.  The  exact  mechanism  through  which  the  helium  thermal 
conductivity  is  affected  by  helium  pressure  is  discussed  later  in  the  thesis. 

The  cooling  water  jacket  is  made  of  Russian  1X18H10T  stainless  steel  and  has  an 
inner  wall  thickness  of  2.5  mm  and  an  outer  wall  thickness  of  1  mm.  A  coiled  stainless 
steel  tube  between  the  jacket  walls  provides  a  spiral  flow  path  for  the  cooling  water.  The 
water  flows  from  the  bottom  to  the  top  of  the  water  jacket  around  the  outside  of  the  coil 
in  the  cooling  "channel"  depicted  in  Figure  4-2.  The  coils  are  spaced  by  35.0  mm  and 
spiral  approximately  6.5  times  from  bottom  to  top.  Cooling  water  temperature  is 
measured  on  the  inlet  and  outlet  piping  exterior  to  the  vacuum  chamber  using  surface 
thermocouples.  Cooling  water  flow  rate  is  measured  in  the  inlet  piping  using  a 
turbine-type  flowmeter  located  exterior  to  the  vacuum  chamber. 
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V.  TEST  STAND  THERMAL  ANALYSIS 


As  discussed  in  Chapter  III,  the  thermionic  cycle  is  a  relatively  inefficient 
conversion  process.  The  TOPAZ-II  TFE  has  a  maximum  efficiency  of  5-6%,  therefore,  at 
least  94%  of  the  heat  generated  by  the  TIS  A  heater  must  be  dissipated  as  a  loss.  Detailed 
thermal  analysis  of  the  test  stand  can  provide  a  better  understanding  of  how  testing  a 
single  TFE  in  the  test  stand  correlates  to  operating  a  TFE  in  the  TOPAZ-II  reactor  core. 
The  developed  thermal  model  also  yields  valuable  information  about  the  limitations  of  the 
test  stand  for  future  testing  and  the  validity  of  past  and  future  test  results. 

A.         THERMAL  MODEL  INPUTS 

The  amount  of  input  power  to  the  test  stand  is  determined  by  measuring  the  input 
current  to  the  TISA  heater  through  an  external  shunt  and  by  measuring  the  voltage  drop 
across  the  input  terminals  of  the  TISA  heater.  Multiplying  the  current  in  amps  by  the 
voltage  in  volts  gives  the  input  power  to  the  TISA  heater  in  watts  and  is  designated  as 
PTISA.  This  TISA  input  power  must  be  corrected  to  calculate  the  amount  of  power  actually 
delivered  to  the  emitter  of  the  TFE.  Correction  is  required  to  take  into  account  the  lead 
and  line  losses  associated  with  the  heater.  Although  Russian  documentation  for 
calculating  the  correction  factor  was  not  available,  the  Russian  technicians  at  TIP  use  a 
correction  factor  of  88%.  Correcting  the  input  power  to  the  TISA  heater  by  this  factor 
provides  the  amount  of  actual  power  delivered  to  the  "active  zone"  of  the  TFE.  This 
active  zone  power  is  designated  as  P^. 

The  output  power  produced  by  the  TFE  is  determined  by  measuring  the  current 
through  the  external  load  bank  and  by  measuring  the  voltage  drop  across  the  upper  and 
lower  emitter  or  collector  terminal  connections  (see  Figure  4-1).  Since  the  voltage  drop  is 
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measured  directly  across  the  interior  terminal  connections  of  either  the  emitter  or 
collector,  no  correction  factor  is  required  to  determine  the  output  voltage  of  the  TFE. 
Multiplying  the  voltage  in  volts  by  the  current  in  amps  gives  the  TFE  output  power  in 
watts  and  is  designated  as  Pout. 

The  only  temperatures  that  are  actually  measured  internally  in  the  test  stand  are  the 
collector  sleeve  temperatures.  These  measurements  occur  1.8  mm  into  the  stainless  steel 
collector  sleeve  from  its  inner  cylindrical  surface  (see  Figure  4-4).  As  discussed  in 
Chapter  IV,  the  thermocouples  in  the  collector  sleeve  provide  axial  temperature 
measurements  along  the  length  of  the  working  section.  There  are  a  total  of  12 
thermocouples  in  the  collector  sleeve  which  are  designated  T56  through  T67.  During  the 
data  collection  runs,  thermocouple  T64  was  inoperative  and  was  disregarded  in  the 
analysis.  Thermocouples  T59  and  T60  and  thermocouples  T61  and  T62  measure  the 
temperature  at  the  same  location  along  the  collector  sleeve.  Relative  locations  of  the 
thermocouples  with  respect  to  the  centerline  of  the  working  section  are  shown  in  Figure 
5-1.  The  lengths  of  the  TISA  heater,  the  thermionic  working  section,  and  the  water  jacket 
and  the  approximate  locations  of  the  cooling  water  inlet  and  outlet  are  also  shown  in 
Figure  5-1. 

In  a  one-dimensional  model,  it  is  necessary  to  take  an  average  of  the  1 1  operating 
collector  sleeve  thermocouples  to  obtain  a  value  to  use  as  an  input  temperature.  Although 
many  variables  effect  the  magnitude  of  the  collector  sleeve  axial  temperature  profile,  the 
profile  shape  is  generally  the  same  under  all  conditions.  Figure  5-2  shows  a  typical  axial 
temperature  profile  at  a  TISA  heater  input  power  of  2400  W.  The  temperature  data 
points  are  shown  as  diamonds  with  T56  at  the  -205  mm  position  and  T67  at  the  205"  mm 
position.  The  dashed  line  represents  a  third-order  polynomial  curve  fit  of  the  data. 
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Figure  5-2.  Typical  Axial  Temperature  Distribution  along  the  Collector  Sleeve. 

Although  the  TISA  heater  provides  a  fairly  uniform  temperature  profile  into  the  emitter, 
several  reasons  inherent  to  the  test  stand  design  cause  the  distribution  to  take  the  shape  as 
shown  in  Figure  5-2.  An  overall  downward  slope  in  the  data  is  expected  due  to  the 
cooling  water  flowing  from  the  bottom  to  the  top  of  the  working  section;  however  the  end 
effects  of  the  upper  and  lower  helium  chambers  cause  the  data  to  not  be  linear  at  the  two 
ends  of  the  graph.  A  straight  average  of  the  1 1  data  points  would  bias  the  temperature 
input  to  the  thermal  model  low  due  mainly  to  the  end  effects.  This  biasing  can  be 
minimized  if  the  end  data  points  (T56  and  T67)  are  not  used  and  an  average  is  taken  only 
of  the  remaining  nine  data  points.  For  example,  the  average  of  the  1 1  data  points  in 
Figure  5-2  is  379.0  K.  Excluding  the  end  data  points  gives  an  average  of  392. 1  K  which  is 
approximately  3.5%  higher.  For  all  calculations  performed  in  this  analysis,  the  average 
collector  sleeve  temperature  was  determined  by  excluding  T56  and  T67.  A  much  more 
detailed  two-dimensional  model  would  be  required  to  understand  the  error  associated  with 
this  approximation. 

32 


The  cooling  water  system  provides  the  last  of  the  measured  test  stand  quantities 
necessary  for  the  thermal  analysis.  Instrumentation  of  the  cooling  water  system  was 
required  in  support  of  this  thesis  since  the  information  on  the  cooling  water  system  was 
not  previously  deemed  necessary  Surface  thermocouples  were  placed  on  the  cooling 
water  inlet  and  outlet  lines,  exterior  to  the  vacuum  chamber  Approximately  one  meter  of 
piping  separates  the  actual  thermocouple  location  from  the  inlet  and  outlet  piping 
penetrations  on  the  test  stand.  No  corrected  factor  was  used  on  these  temperature 
measurements  for  any  potential  piping  losses  between  the  test  section  and  the  exterior  of 
the  vacuum  chamber  because  the  only  method  of  heat  transfer  from  the  water  piping 
would  be  by  radiation  heat  transfer  to  the  interior  of  the  vacuum  chamber.  Since  the 
maximum  water  temperature  inside  the  piping  is  less  than  40  °C  over  all  the  data  runs, 
radiation  heat  transfer  is  not  significant.  The  water  piping  is  also  in  thermal  contact  with 
the  vacuum  chamber  which  provides  another  potential  heat  loss  by  conduction  heat 
transfer  to  the  vacuum  chamber  itself.  Some  uncertainty  could  have  been  removed  if  the 
thermocouples  had  been  placed  on  the  piping  inside  the  vacuum  chamber  at  the  entrance 
and  exit  of  the  test  section.  This  was  not  done  for  any  of  the  data  runs  used  in  this 
analysis  due  to  time  constraints  and  the  complexity  of  removing  the  vacuum  chamber.  An 
average  bulk  water  temperature  in  the  water  jacket  can  be  estimated  by  taking  an  average 
of  the  inlet  and  outlet  water  temperatures. 

A  turbine-type  flow  meter  was  also  installed  on  the  inlet  cooling  water  line  prior  to 

the  data  runs.   The  accuracy  of  the  flow  meter  was  physically  verified  to  be  within  2%  of 

the  indicated  flow  rate  over  the  range  used  during  the  data  runs.    Initially  the  data  runs 

were  taken  with  the  cooling  water  isolation  values  fully  open  which  was  the  normal 

operating  procedures  for  the  test  stand.   However,  with  the  valves  fully  open,  the  change 

in  cooling  water  temperature  between  the  inlet  and  the  outlet  was  less  than  5  °C.  Since 
the  accuracy  of  the  surface  thermocouples  is  only  ±  1  °C,  a  larger  temperature  change  is 
desired  to  reduce  the  percent  error.  Around  halfway  through  the  first  set  of  data  runs,  the 

flow  rate  was  throttled  down  to  approximately  0.7  gpm.  This  gave  a  temperature  change 
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between  10  to  15  °C  over  most  of  the  data  runs.    All  of  the  data  collected  during  the 
second  set  of  data  runs  were  collected  at  a  flow  rate  of  approximately  0.7  gpm. 

B.         THERMAL  MODEL 

Difficulty  was  encountered  during  the  thermal  analysis  of  the  test  stand  due  to  the 
small  number  of  input  parameters  and  the  large  number  of  "unknowns."  In  order  to  better 
understand  the  scope  of  the  problem,  a  thermal  analysis  flowchart  was  developed.  Figure 
5-3  shows  the  test  stand  thermal  analysis  flowchart.  Data  collected  from  the  test  stand  are 
shown  as  circle  and  are  used  as  inputs  into  the  thermal  model.  From  this  flowchart  it  can 
be  seen  that  the  problem  divides  itself  into  the  following  three  sections:  (1)  thermal 
network  calculations;  (2)  water  channel  heat  removal  calculations;  and  (3)  coupling  the 
thermal  network  calculations  with  the  TOPAZ-II  Thermionic  Transient  Analysis  Model 
(TIT  AM).  TIT  AM  is  a  computer  model  developed  by  the  University  of  New  Mexico  to 
simulate  the  operation  of  the  entire  TOPAZ-II  reactor  system  under  steady-state  and 
transient  conditions  (Paramonov  and  El-Genk,  1994).  The  basis  for  the  entire  analysis  is 
the  thermal  network. 

1.  Thermal  Network 

The  one-dimensional  losses  through  the  test  stand  are  modeled  using  a  network 
method.  Figure  5-4  shows  a  simplified  representation  of  the  resistance  network  used  for 
thermal  modeling  the  test  stand.  The  network  starts  at  the  outside  wall  of  the  cylindrical 
emitter  and  runs  outward  to  the  centerline  of  the  water  channel.  The  three  major  modes 
of  heat  transfer  encountered  throughout  the  network  are  radiation,  conduction,  and  forced 
convection.  Radiation  heat  transfer  was  considered  a  heat  transfer  mode  through  all  of  the 
test  stand  gaps  when  the  thermal  network  was  initially  developed.  It  is  shown  later  in  this 
section  that  radiation  heat  transfer  plays  an  insignificant  role  in  the  heat  transferred 
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through  the  unregulated  and  regulated  helium  gaps,  however,  radiation  heat  transfer  does 
play  an  important  role  in  the  cesium  interelectrode  gap.  Conduction  heat  transfer  occurs 
through  all  of  the  metal  components  and  through  the  gaps  in  the  test  stand;  including  the 
cesium  interelectrode  gap,  the  two  unregulated  helium  gaps,  and  the  regulated  helium  gap. 
Forced  convection  occurs  into  the  cooling  water  flowing  through  the  water  channel. 

The  basic  heat  transfer  equation  used  for  conduction  heat  transfer  through  a 
cylindrical  wall  is  found  in  Wong  (1977).  Equation  5.1  gives  the  equation  used  for  the 
metal  components  of  the  test  stand: 


Q  =  ^-  (Thot  -  Tcold)  (5.1) 

Kc 


where:  Rc  =  thermal  resistance  (K/W) 
Thot  =  hot  wall  temperature  (K) 
Tcoid =  c°ld  wall  temperature  (K) 


The  thermal  resistance,  Rc,  is  calculated  for  hollow  cylinders  using  Equation  5.2: 


where:  r0  =  outer  wall  radius  (m) 
rt  =  inner  wall  radius  (m) 
z  =  cylinder  length  (m) 
k  =  material  thermal  conductivity  (W/mK) 


For  the  test  stand,  Thot  will  be  the  temperature  at  the  inner  cylindrical  surface  and  Tcold  will 
be  the  temperature  at  the  outer  cylindrical  surface  of  the  component  of  interest.  For  the 
thermal  model  calculations,  the  temperature  subscript  will  indicate  the  node  at  which  the 
temperature  is  measured. 
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A  cylindrical  length  of  30.0  cm  is  used  in  all  of  the  calculations.     This  length 

corresponds  to  the  total  length  of  the  TISA  heater;  however,  as  can  be  seen  in  Figure  5-1, 

the  length  of  the  TISA  heater,  the  length  of  the  water  jacket,  and  the  length  of  the  TFE 

working  section  are  all  different.   Since  there  is  a  20%  increase  in  cylinder  length  from  the 

TISA  heater  length  to  the  TFE  working  section  length,  the  thermal  resistance  also 
increases  by  20%.   Although  this  only  results  in  a  few  degree  Kelvin  difference  in  the  AT 
across  a  single  metal  component,  the  cumulative  effect  of  this  difference  can  lead  to  large 

calculational  uncertainty.     The  shortest  of  the  possible  cylinder  lengths  is  used  in  the 

analysis  since  it  gives  a  more  conservative  result. 

The  thermal  conductivities  of  all  the  different  test  stand  components  are  graphed  as 

a  function  of  temperature  in  Appendix  A.    Since  the  thermal  conductivity  is  a  function  of 

temperature  and  the  metal  components  have  a  decreasing  temperature  distribution  from 

the  inner  wall  to  the  outer  wall,  an  estimate  of  the  thermal  conductivity  for  the  entire 

thickness  must  be  made.     When  entering  the  thermal  model,  only  one  temperature  is 

known;   therefore,   it  is  not  possible  to  determine  the  average  temperature  of  the 

component.    For  the  first  iteration  through  the  network,  the  thermal  conductivity  of  the 

component  is  calculated  based  on  the  one  known  temperature.  Because  metal  is  a  very 
good  conductor  of  heat,  the  thermal  conductivity  is  fairly  constant  over  the  small  AT 
effective  across  the  metal.    Based  on  calculations,  using  the  one  known  temperature  to 

estimate  the  thermal  conductivity  results  in  only  a  0.4%  difference  than  if  the  thermal 
conductivity  were  recalculated  with  the  lower  temperature  found  after  calculating  the  AT. 
Therefore,  less  than  a  three  degree  Kelvin  error  is  induced  into  the  network  and  this  makes 

recalculation  of  the  thermal  conductivity  unnecessary. 

The  basic  equation  used  for  heat  transfer  for  the  gas  filled  gaps  is  found  in  El-Oenk 
et  al.  (1993).  This  equation  accounts  for  both  radiation  and  conduction  heat  transfer 
through  the  gap  and  is  given  in  Equation  5.3: 
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Q  =  zr~  (Th0,  -  Tcold  (5.3) 

Kr/c 


where:  R,  c  =  radiation/conduction  thermal  resistance  (K/W) 


The  radiation/conduction  thermal  resistance,  R^,  is  calculated  using  Equation  5.4: 

ln(r)      k 

R'/c  =  i     <  ,{T  +  CT£eg<Thot+Tcoid)(Thot+TCO|d)}-1  (5.4) 

27iz(r0-rj)    d 

where:  kg  =  gas  thermal  conductivity  (W/mK) 
d  =  gap  width  (m)  =  (ro  -  r,) 

a  =  Stefan-Boltzmann  constant  =  5.67xl0"8  W/m:K4 
eeff  =  effective  emissivity 


The  thermal  conductivity  for  helium  as  a  function  of  temperature  is  graphed  in  Appendix 
A.  According  to  Reid  et  al.  (1977),  the  thermal  conductivity  of  all  gases  can  also  be 
affected  by  the  pressure  in  the  gap,  although  for  pressures  greater  than  approximately  1 
torr,  the  thermal  conductivity  only  increases  by  about  1%  for  a  pressure  change  of  760 
ton.  Therefore,  for  the  unregulated  helium  gap  pressures  of  200-300  torr,  the  thermal 
conductivity  is  considered  a  function  of  temperature  only.  This  estimation  is  valid  because 
the  mean  free  path  for  collisions  between  atoms  is  much  less  than  the  mean  free  path  for 
collision  between  an  atom  and  the  walls;  therefore,  an  atom  does  not  have  to  travel  far 
before  transferring  its  energy  to  another  atom  and  eventually  reaching  thermal  energies. 
The  pressure  maintained  in  the  regulated  helium  gap  is  around  the  lower  end  of  where  this 
approximation  is  considered  valid.  According  to  the  experimental  data  collected,  the 
thermal  conductivity  of  helium  in  the  1-10  torr  range  is  greatly  affected  by  pressure.  This 
unusual  phenomenon  is  discussed  in  detail  later  in  this  chapter. 

The  effective  emissivity  can  be  calculated  by  assuming  that  the  view  factors  from 
both  surfaces  are  equal  to  1,  and  the  two  surfaces  areas  are  also  equal  by  Equation  5.5 
(Weltyetal.,  1984): 
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8'ff=l J (5.5) 

—  +  —  -  1 

Ej        £o 


where:  6j  =  inner  surface  emissivity 
e0  =  outer  surface  emissivity 


Due  to  the  extremely  small  dimensions  and  gap  size,  both  of  these  assumptions  are  valid. 

The  basic  equation  used  for  forced  convection  heat  transfer  into  the  cooling  water 
is  found  in  Wong  (1977).  The  equation  used  for  calculating  the  heat  transfer  into  the  test 
stand  cooling  water  is  based  on  Newton's  Law  of  Cooling  and  is  given  in  Equation  5.6: 

Q  =  7T~  (Twall  -  Tbulk)  (5  6) 


where:  R^  =  forced  convection  thermal  resistance  (K/W) 
Tw,!!  =  cooling  water  channel  wall  temperature  (K) 
Tbuik =  bulk  cooling  water  temperature  (K) 


The  forced  convection  thermal  resistance,  R^.,  is  calculated  using  Equation  5.7: 

In(£) 

RfC=       :  Fl    rr  (5.7) 

27tz(r0-ri)h 

where:  h  =  convective  heat  transfer  coefficient  (W/m2K) 


The  convective  heat  transfer  coefficient  is  a  complicated  function  of  the  cooling  channel 
geometry,  the  thermophysical  properties  of  the  cooling  water,  and  principally  of  the  flow 
characteristics  through  the  channel.  This  coefficient  and  the  complexities  encountered  in 
its  definition  are  discussed  in  detail  later  in  this  chapter. 
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An  example  of  the  data  set  used  in  the  thermal  network  calculations  is  given  in 
Table  5-1  These  data  are  used  in  the  following  discussion  to  explain  the  thermal  model 
calculational  procedure.  The  data  collected  between  14-18  August,  1994,  in  support  of 
this  thesis  are  tabulated  in  Appendix  B  The  data  presented  in  Table  5-1  were  taken  on 
August  17,  1994.  As  can  be  seen  there,  the  TFE  efficiency  is  less  than  4%.  This  is  mainly 
due  to  the  interelectrode  cesium  pressure  not  being  at  its  optimum  value  for  an  input 
power  of  3000  W.  At  this  power  level,  an  interelectrode  cesium  pressure  of 
approximately  1  torr  would  increase  the  TFE's  efficiency  closer  to  5%.  However,  for 
thermal  modeling  purposes,  optimum  TFE  performance  is  not  important.  Cesium 
optimization  will  be  covered  in  detail  in  the  companion  thesis  written  by  LT  J.  R.  Venable 
that  will  be  available  after  March  of  1995. 


Ptisa  W 

3,412 

Minimum  Collector 
Sleeve  Temp.  (C) 

354.3 

PaJW) 

3,003 

Maximum  Collector 
Sleeve  Temp.  (C) 

515.6 

P-  (W) 

109 

Average  Collector 
Sleeve  Temp.  (C)  T6 

453.7 

Q,«(W) 

2,894 

Cooling  Water  Flow 
rate  (gpm) 

0.69 

Efficiency  (%) 

3.63 

Cooling  Water  Mass 
Flow  rate  (kg/s) 

0.04 

Cs  Pressure 
(torr) 

0.5 

Cooling  Water  Inlet 
Temp.  (C) 

16.6 

Regulated  He 
Gap  Press  (torr) 

10 

Cooling  Water  Outlet 
Temp.  (C) 

32.2/32.1 

Table  5-1.  Example  Data  Set  for  the  Thermal  Model. 
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Since  the  only  temperature  measured  internally  to  the  test  stand  is  the  collector 
sleeve  temperature,  the  thermal  network  is  entered  at  node  T6.  From  this  node, 
calculations  can  be  made  to  determine  the  temperature  distribution  either  inward  toward 
the  center  of  the  TFE  or  outward  toward  the  water  channel.  In  order  to  make  these 
temperature  calculations  using  basic  heat  transfer  equations,  the  amount  of  heat 
transferred  out  through  the  test  stand  must  be  calculated.  This  transferred  heat  is  equal  to 
the  portion  of  the  TISA  input  power  that  is  not  converted  to  electricity  in  the  TFE  and  is 
designated  as  Qloss.  Equation  5.8  shows  the  definition  for  Qloss: 

Qloss  =  Paz  -  Pout  (58) 

In  a  one-dimensional  model  for  the  test  stand,  Q|0M  will  travel  only  radially  outward 
through  the  test  stand  components  to  be  carried  away  by  the  cooling  water  flowing 
through  the  water  jacket. 

Using  the  node  6  temperature  (T6)  and  the  amount  of  heat  lost  by  the  system  (Qloss) 
as  a  starting  point,  the  thermal  model  is  presented  in  two  sections:  (1)  calculations  from 
the  collector  sleeve  to  the  water  channel,  and  (2)  calculations  from  the  collector  sleeve  to 
the  emitter.  Detailed  discussions  about  the  critical  resistances  in  the  network  are  also  be 
presented.  The  calculations  made  with  the  data  from  Table  5-1  are  given  in  Appendix  C. 

a.  From  the  Collector  Sleeve  to  the  Water  Channel 

The  first  resistance  encountered  in  the  thermal  network  is  a  1.2  mm  thick 
stainless  steel  collector  sleeve.  The  thermal  conductivity  of  the  stainless  steel  is 
determined  using  Figure  A-3.  This  thermal  conductivity  is  then  used  to  calculate  the 
thermal  resistance  of  the  stainless  steel  using  Equation  5.2.   Using  the  calculated  Qloss  as 
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the  variable  Q,  Equation  5.1  can  now  be  solved  for  Tcold,  which  corresponds  to  T7,  since  all 
of  the  other  equation  variables  have  been  defined. 

A  0.05  mm  unregulated  helium  gap  is  the  next  resistance  in  the  thermal 
network.  The  thermal  conductivity  of  this  continuum  helium  gas  (gap  pressure  of 
approximately  200  torr)  is  calculated  using  Figure  A-4  The  effective  emissivity  is  then 
calculated  using  Equation  5.5.  This  equation  requires  the  stainless  steel  emissivity  which 
is  found  from  Figure  A-6  and  the  copper  emissivity  which  is  found  from  Figure  A-7.  Both 
of  these  emissivities  are  fairly  constant  over  the  temperature  range  of  interest.  Using 
Equation  5.4  in  combination  with  Equation  5.3,  the  value  for  T8  can  be  solved  for 
numerically.  If  the  contribution  of  the  radiation  heat  transfer  is  ignored  through  the  helium 
gap,  and  T8  is  solved  for  based  on  conduction  heat  transfer  through  the  gap  only,  the 
answers  for  Tg  are  the  same  within  four  significant  digits.  For  this  reason,  radiation  heat 
transfer  through  the  unregulated  and  regulated  helium  gaps  can  be  disregarded  and  the 
resistance  network  can  be  solved  based  on  conduction  heat  transfer  only. 

The  next  resistance  in  the  network  is  the  copper  conductor.  Again 
Equations  5.2  and  5.1  are  solved  to  determine  the  temperature  on  the  outer  cylindrical 
surface  of  the  copper  conductor,  which  corresponds  to  T9.  The  regulated  helium  gap  is 
the  next  resistance  in  the  thermal  network. 

b.  Regulated  Helium  Gap 

As  discussed  earlier  in  this  chapter,  the  pressure  of  the  helium  in  the 
regulated  gap  is  maintained  between  1-10  torr,  depending  on  the  input  power  level,  to 
hold  the  collector  sleeve  temperature  fairly  constant.  Since  the  collector  sleeve 
temperatures  are  measured  at  approximately  the  same  radial  location  that  the  liquid  metal 
cooling  flows  through  the  actual  TOPAZ-II  core,  holding  this  temperature  constant  better 
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simulates  the  actual  conditions  in  the  core.  The  temperature  at  this  radial  position  in  the 
TOPAZ-II  core  is  relatively  constant  due  to  the  high  heat  capacity  of  the  liquid  metal 
cooling  and  the  high  heat  rejection  capability  of  the  radiator. 

Due  to  the  low  gap  pressure,  the  helium  is  rarefied  which  indicates  that  the 
mean  free  path  for  collisions  between  atoms  is  much  larger  than  the  mean  free  path  for 
atom  to  wall  collisions.  Paul  Agnew  of  AEA  Technology,  United  Kingdom,  currently 
assigned  to  the  TOPAZ-II  Project  and  working  in  the  area  of  materials  science  and 
thermionics,  has  completed  a  preliminary  analysis  on  the  pressure  dependence  of  helium 
thermal  conductivity  at  these  low  pressures.  The  initial  results  of  his  work  indicated  that 
helium  thermal  conductivity  should  approach  its  continuum  value  at  a  gap  pressure  of 
approximately  5  torr  vice  the  1  torr  discussed  in  Reid  et  al.  (1977).  However, 
experimental  results  indicate  that  the  helium  thermal  conductivity  should  be  around  30% 
of  its  continuum  value  for  a  gap  pressure  of  10  torr.  Detailed  Monte  Carlo  computer 
programming  and  investigation  of  gas  atom  scattering  mechanisms  will  be  required  to 
determine  a  more  accurate  model  of  the  regulated  gap.  (Agnew,  1994) 

Although  a  more  exact  solution  is  not  possible  without  extensive  work,  the 
regulated  gap  helium  thermal  conductivity  can  be  better  approximated  using  the  empirical 
data  collected  on  the  test  stand.  Since  the  regulated  helium  gap  is  separated  from  the 
cooling  water  channel  by  only  2.0  mm  of  stainless  steel,  the  temperature  of  the  inner  wall 
of  the  cooling  channel  is  directly  dependent  on  the  regulated  helium  gap  thermal 
conductivity.  Using  the  helium  continuum  thermal  conductivity  value  calculated  to  be 
0.276  W/mK  for  this  example,  gives  a  cooling  channel  inner  wall  temperature  of  around 
325  °C.  The  cooling  water  is  pressurized  on  average  to  a  pressure  not  greater  than  70 
psig.  Since  the  cooling  system  provides  water  to  many  different  components,  the  actual 
pressure  is  dependent  on  the  cooling  system  load  conditions.  This  cooling  water  pressure 
gives  a  initiation  boiling  point  temperature  of  approximately  140  °C.  The  wall  temperature 
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calculated  using  the  helium  continuum  thermal  conductivity  is  around  185  °C  greater  than 

this  boiling  point  temperature.    Based  on  an  cooling  water  outlet  temperature  of  32  °C,  it 

is  assumed  that  the  water  does  not  boil  in  the  cooling  channel.   Figure  5-5  shows  a  graph 

of  the  cooling  water  channel  inner  wall  temperature  dependence  on  the  regulated  helium 
gap  thermal  conductivity.  If  the  wall  temperature  is  bracketed  around  80  ±  30  °C,  the 
regulated  helium  gap  thermal  conductivity  ranges  between  0.073-0.087  W/mK.     This 

range  is  around  30%  of  the  calculated  continuum  value. 

This  calculated  range  of  helium  gap  thermal  conductivities  can  be  used 
along  with  Equation  5. 1  to  calculated  a  range  of  node  10  temperatures.  Using  these  node 
10  temperatures,  the  2.0  mm  thick  stainless  steel  water  jacket  thermal  conductivity  can  be 
determined  using  Figure  A-3.    The  temperature  at  the  inner  wall  of  the  cooling  water 
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Figure  5-5.  Cooling  Water  Channel  Inner  Wall  Temperature  Dependence  on 
the  Regulated  Helium  Gap  Thermal  Conductivity. 
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channel,  corresponding  to  the  node  1 1  temperature,  is  then  calculated  using  Equations  5.2 
and  5.1.   These  calculated  wall  temperatures  range  from  64-122  °C. 

Using  an  estimated  convective  heat  transfer  coefficient  and  an  estimated 
heat  transfer  area  (the  method  used  to  estimate  these  parameters  is  discussed  later  in  this 
chapter),  the  bulk  water  temperature  can  be  estimated  using  Equations  5.7  and  5.6.  These 
calculations  yield  a  range  of  bulk  water  temperatures  of  6-92  °C.  Using  an  average  of  the 
measured  inlet  and  outlet  cooling  water  temperatures  gives  a  bulk  temperature  of 
approximately  24  °C. 

c.  From  the  Collector  Sleeve  to  the  Emitter 

The  calculations  for  the  inward  portion  of  the  thermal  network  are  similar 
to  the  calculations  made  for  the  outward  portion.  The  first  resistance  encountered  is  the 
inner  portion  of  the  stainless  steel  collector  sleeve.  The  same  thermal  conductivity  value  is 
used  for  the  collector  sleeve  as  was  used  in  the  outward  portion.  Using  Equations  5.2  and 
5. 1,  the  node  5  temperature  can  be  found.  This  corresponds  to  T5. 

The  next  resistance  in  the  network  is  another  0.05  mm  unregulated  helium 
gap.  Radiation  heat  transfer  is  ignored  and  the  heat  transfer  calculations  using  Equations 
5.2  and  5.1  provide  the  temperature  at  node  4.  Node  4  corresponds  to  the  outer  wall  of 
the  layer  of  alumina  insulation  on  the  collector.  The  thermal  conductivity  of  alumina  is 
found  using  Figure  A- land  the  same  conduction  heat  transfer  calculations  are  completed 
to  calculate  the  node  3  temperature,  T3.  This  temperature  is  the  interface  temperature 
between  the  collector  and  the  alumina  insulation. 

The  temperature  at  node  2,  T2,  can  be  similarly  calculated  using  the 
thermal  conductivity  of  polycrystalline  molybdenum  found  from  Figure  A-5.      This 
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temperature  corresponds  to  the  inner  surface  of  the  collector  and  is  an  extremely 
important  measurement  to  determine  thermionic  performance  The  cesium  interelectrode 
gap  is  the  next  resistance  in  the  network. 

d.   Cesium  Interelectrode  Gap 

The  cesium  interelctrode  gap  heat  transfer  mechanism  is  another  critical 
resistance  in  the  network  that  cannot  be  easily  calculated  with  only  basic  heat  transfer 
equations.  Conduction  and  radiation  heat  transfer  both  play  an  important  part  in  the 
overall  heat  transfer  through  the  interelectrode  gap.  The  thermal  conductivity  of  the 
cesium  in  the  gap  depends  on  the  gases  vapor  pressure,  the  emitter  surface  temperature, 
the  collector  surface  temperature,  and  the  gap  size.  (Paramonov  and  El-Genk,  1994) 

Radiation  heat  transfer  cannot  be  ignored  through  the  interelectrode  gap 
since  the  temperature  of  the  emitter  is  high  enough  to  make  this  mode  of  heat  transfer 
significant  (on  the  order  of  1800  K).  The  effective  emissivity  for  the  gap  cannot  be 
calculated  as  presented  for  the  helium  gaps,  but  must  be  determined  empirically.  A 
detailed  discussion  on  thermal  modeling  of  the  cesium  interelectrode  gap  and  coupling  the 
thermal  model  with  a  thermionic  emission  model  can  be  found  in  Paramonov  and  El-Genk 
(1994). 

2.  Water  Channel  Calculations 

In  order  to  calculate  the  thermal  resistance  using  Equation  5.7,  the  convective  heat 
transfer  coefficient  is  required.  Calculating  the  heat  transfer  coefficient  for  the  Water 
channel  is  complicated  by  the  flow  channel  geometry  and  by  the  low  water  flow  rates. 
Figure  5-6  shows  the  cooling  water  jacket  details  including  the  flow  guidance  coil.  This 
flow  guidance  coil  is  a  stainless  steel  tube  with  an  outer  diameter  of  10.0  mm  that  fits 
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Flow    Guidance    Tube    in    the    Water    Jacket 
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Figure  5-6.  Cooling  Water  Jacket  Details. 
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inside  of  the  water  jacket  to  direct  the  flow  in  a  spiral  path  to  facilitate  forced  convection 
heat  transfer.  The  cooling  water  flows  around  the  outside  of  the  flow  guidance  coil  with 
no  flow  through  the  coil  itself. 

As  shown  in  Figure  5-6,  the  flow  channel  cross-sectional  area  is  calculated  by 
determining  the  cross-sectional  area  of  a  rectangle  with  dimensions  of  35x10  mm,  and 
subtracting  off  the  cross-sectional  area  of  the  flow  guidance  coil  (diameter  of  10  mm). 
This  gives  a  cross-sectional  channel  flow  area  of  2. 7x10^  m:.  The  equivalent  hydraulic 
diameter  is  calculated  using  the  definition  from  Wong  (1977),  to  be  0.01 1  m. 

Following  the  definitions  and  equations  presented  by  El-Genk  and  Rao  (1986),  the 
Reynolds  number  is  calculated  using  this  geometry  and  the  given  example  flow  rate  to  be 
1480  This  Reynolds  number  indicates  that  the  flow  falls  into  the  combined  laminar  flow 
regime.  The  convective  heat  transfer  coefficient  can  be  defined  from  El-Genk  and  Rao 
(1986)  and  is  given  in  Equation  5.9: 

.  _Nuk 

h_~DT  <59> 

where:  Nu  =  Nusselt  number 

k  =  water  thermal  conductivity  (W/mK) 
De  =  equivalent  diameter  (m) 

In  the  combined  laminar  regime,  the  Nusselt  number  is  dependent  on  the  thermophysical 
properties  (defined  by  the  Graetz  number)  and  the  buoyant  forces  (defined  by  the  Rayleigh 
number).  This  flow  regime  is  immensely  difficult  to  analyze  with  existing  correlations 

Although  the  convective  heat  transfer  coefficient  cannot  be  determined  exactly, 
several  existing  forced  convection  correlations  can  provide  a  range  of  possible 
coefficients.    The  upper  limit  is  determined  using  the  Dittus-Boelter  correlation  for  fully 
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developed  turbulent  flow  in  a  channel.  This  gives  a  convective  heat  transfer  coefficient  of 
1012  W/m:K.  The  lower  limit  is  found  by  the  Seider-Tate  correlation  for  developing 
laminar  flow  This  correlation  gives  a  convective  heat  transfer  coefficient  of  528  VV/nrK. 
Therefore,  the  actual  coefficient  will  fall  between  the  range  of  528-1012  W/m:K  and  these 
limits  can  then  be  used  in  the  resistance  network  to  estimate  the  drop  between  the  cooling 
channel  wall  temperature  and  the  bulk  water  temperature. 

Using  the  cooling  water  inlet  and  outlet  temperatures  along  with  the  measured 
cooling  water  flow  rate,  the  amount  of  heat  carried  away  by  the  water  can  be  calculated 
using  Equation  5.10: 

Q  =  rh  cp  (Tout  -  Tin)  (5.10) 


where:  rh  =  cooling  water  mass  flow  rate  (kg/s) 
cp  =  specific  heat  capacity  of  water 
Tout  =  cooling  water  outlet  temperature  (K) 
Tm  =  cooling  water  inlet  temperature  (K) 


The  heat  carried  away  by  the  water,  Q,  should  be  equal  to  the  heat  lost  by  the  TEE,  Qloss,  if 
all  of  the  excess  heat  traveled  one-dimensionally  out  of  the  test  section  to  be  carried  away 
by  the  cooling  water.  For  the  example  data,  the  heat  carried  away  by  the  cooling  water  is 
approximately  2840  W.  This  is  less  than  2%  different  from  the  2894  W  lost  by  the  TFE. 
As  can  be  seen  from  the  experimental  data  presented  in  Appendix  B,  most  of  the  data 
collected  for  the  heat  removed  by  the  cooling  water  falls  within  10%  of  the  heat  lost  by 
the  TFE.  The  percent  error  tends  to  increase  as  the  input  power  level  decreases.  This  is 
an  expected  result  since  the  difference  between  the  inlet  and  outlet  cooling  water 
temperature  decreases  as  the  power  level  decreases  causing  an  increase  in  the  expected 
error. 
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By  setting  Equation  5  10  equal  to  Equation  5.6,  the  convective  heat  transfer 
coefficient  can  be  calculated  based  on  a  heat  transfer  surface  area  Once  again  difficulty 
arises  in  determining  the  exact  heat  transfer  surface  area.  From  Figure  5-6  it  can  be  seen 
that  even  though  no  cooling  water  flows  through  the  flow  guidance  coil,  the  stainless  steel 
coil  is  in  thermal  contact  with  the  stainless  steel  water  jacket  The  maximum  possible  heat 
transfer  surface  area  can  be  determined  by  summing  the  inner  water  jacket  cylindrical 
surface  area,  the  outer  water  jacket  cylindrical  surface  area,  and  the  total  surface  area  of 
the  flow  guidance  coil.  The  minimum  possible  heat  transfer  surface  area  is  equal  to  the 
cylindrical  surface  area  of  the  inner  water  jacket.  These  calculations  are  shown  in 
Appendix  D  and  give  a  range  of  heat  transfer  area  between  0.067-0.203  m:.  Figure  5-7 
shows  a  graph  of  the  heat  transfer  coefficient  plotted  against  the  heat  transfer  area.    Since 
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Figure  5-7.  Cooling  Water  Heat  Transfer  Coefficient  Dependence  on  the 
Heat  Transfer  Area. 
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the  cooling  water  bulk  temperature  is  known,  the  cooling  channel  wall  temperature  is 
varied  from  45  °C  to  95  °C.  It  can  be  seen  from  the  figure  that  a  wall  temperature  of 
greater  than  95  °C  would  fall  outside  of  the  heat  transfer  coefficient/heat  transfer  area 
bands.  This  information  can  be  used  along  with  Figure  5-5  to  further  limit  the  regulated 
helium  gap  thermal  conductivity  used  in  the  thermal  resistance  network. 

3.  Thermal  Network  Coupling  With  TITAM 

The  TITAM  computer  code  provides  a  means  of  simulating  the  operation  of  a 
TOPAZ-II  TFE  in  the  reactor  core  under  both  steady-state  and  transient  conditions 
(Paramonov  and  El-Genk,  1994).  Although  the  test  stand  is  constructed  differently  than 
the  actual  core  and  the  thermal  and  electrical  coupling  of  the  TFE  with  the  actual  reactor 
systems  does  not  exist,  the  location  of  the  collector  sleeve  thermocouples  in  the  test  stand 
does  provide  a  means  of  entering  TITAM  and  having  TITAM  determine  the  operating 
characteristics  of  the  TFE. 

The  required  inputs  into  the  electrically  heated  TFE  component  operation  portion 
of  TITAM  are  thermal  input  power,  inlet  NaK  coolant  temperature,  and  NaK  coolant  flow 
rate.  The  thermal  input  power  to  the  TFE  is  equal  to  the  active  zone  power,  Ptt,  used  in 
the  thermal  network.  The  average  collector  sleeve  temperature,  T6,  may  be  used  as  the 
NaK  coolant  inlet  temperature,  and  a  typical  NaK  flow  rate  of  1.3  kg/s  is  used  as  the  last 
required  input  to  TITAM.  From  the  inputted  data,  TITAM  outputs  the  radial  temperature 
distribution  through  a  TFE  in  the  reactor  core  and  provides  the  output  operating 
characteristics.  The  TITAM  program  was  executed  using  the  inputs  from  Table  5-1.  The 
steady-state  output  and  the  relevant  data  from  the  thermal  analysis  are  presented  in  Table 
5-2. 
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Parameter 

TITAM  Output 

Thermal  Analysis 
Results 

Thermal  Power  to 
TFE  (W) 

3,003 

3,003 

TFE  Load  Voltage 
(V) 

0.81 

1.38 

Load  Current 
(A) 

18438 

78.89 

Load  Electric  Power 
(W) 

149.58 

108.87 

TFE  Efficeincy 
(%) 

4.98 

3.63 

Cesium  Reservior 
Temperature  (K) 

576 

528 

Emitter  Surface 
Temperature  (K) 

1,835.57 

- 

Collector  Surface 
Temperature  (K) 

765.71 

763.93 

Insulator 
Temperature  (K) 

763.75 

761.88 

Coolant  Average 
Temperature  (K) 

727.58 

726.83 

Table  5-2.  Comparison  of  TITAM  output  to  Thermal  Analysis  Results 

The  discrepancy  in  the  TFE  output  power  could  be  the  result  of  the  differences 
between  operating  a  TFE  in  the  test  stand  and  simulating  the  operation  of  a  TFE  in  the 
reactor  core.  Another  possible  source  of  error  are  the  initial  and  operating  conditions  that 
TITAM  uses  in  its  calculations  being  different  from  the  actual  test  stand  conditions. 
TITAM  has  a  user  friendly  interface  through  two  input  files  which  can  be  changed  to 
reflect  various  initial  and  operating  conditions  in  the  reactor  core.  This  includes  the  ability 
to  vary  the  steady-state  interelectrode  gap  cesium  pressure. 
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TITAM  initially  uses  a  steady-state  cesium  pressure  of  approximately  2.0  torr 
(based  on  a  cesium  reservoir  temperature  of  576  K).  The  input  files  were  not  altered  prior 
to  running  the  example  data  set  through  TITAM  to  obtain  the  above  results.  When 
TITAMs  initial  steady-state  cesium  pressure  was  changed  to  reflect  the  actual  operating 
pressure  in  the  test  stand  of  0.5  torr  (based  on  a  cesium  reservoir  temperature  of  528  K), 
the  TFE  output  power  increased  from  151  W  to  188  W.  Although  TITAM  is  being  used 
for  a  purpose  that  it  was  not  designed  for,  this  trend  is  contrary  to  the  experimental  data 
collected. 

The  TITAM  temperature  distribution  correlates  well  with  the  thermal  analysis 
mainly  due  to  these  calculations  being  based  on  simple  heat  transfer  principles.  This 
indicates  that  TITAM  could  be  used  to  estimate  the  emitter  surface  temperature  based  on 
test  stand  thermal  analysis  data. 

In  order  to  understand  the  discrepancies  encountered  by  trying  to  couple  the 
experimental  test  stand  data  with  TITAM,  a  detailed  examination  of  how  TITAM  actually 
performs  it  calculations  and  the  assumptions  under  which  it  operates  is  required. 
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VI.  CONCLUSIONS 


The  TOPAZ-II  TFE  test  stand  provides  an  effective  means  to  investigate  the 
performance  of  an  individual  TFE  without  the  interfering  effects  and  complexity  of  the 
integrated  TOPAZ-II  system.  A  detailed  thermal  analysis  is  completed  to  better  predict 
the  behavior  of  a  single-cell  TFE  and  to  understand  the  thermal  losses  associated  with  the 
test  stand  components. 

Prior  to  developing  the  thermal  model,  extensive  research  of  the  test  stand 
construction  and  materials  was  necessary.  Due  to  the  large  volume  of  Russian 
documentation  and  the  difficulty  in  translating  highly  technical  data  and  engineering 
drawings,  the  development  of  test  stand  cross-sectional  drawings  and  the  determination  of 
component  materials  and  material  properties  were  crucial  to  the  thermal  analysis. 

A  thermal  resistance  network  has  been  utilized  to  determine  the  thermal  loss 
mechanisms  through  each  of  the  test  stand  components.  Conductive  heat  transfer  through 
the  metal  components  accounts  for  a  very  small  portion  of  the  overall  temperature  drop 
from  the  emitter  to  the  water  channel.  For  the  data  set  used  in  the  example  calculations, 
only  7%  of  the  estimated  temperature  drop  from  the  inner  surface  of  the  collector  sleeve 
to  the  bulk  water  temperature  occurred  in  the  metal  components. 

Both  radiation  and  conduction  heat  transfer  occur  through  the  regulated  and 

unregulated  helium  gaps,  but  calculations  indicate  that  radiation  heat  transfer  does  not 

contribute  significantly  to  the  total  heat  transferred  through  these  gaps.    The  pressure  in 

the  regulated  helium  gap  is  varied  during  power  changes  to  maintain  the  collector  sleeve 
temperature  relatively  constant.  Due  to  the  very  low  pressure  maintained  in  this  gap  (<  10 
torr),  the  helium  is  rarefied  and  the  thermal  conductivity  of  the  gas  is  not  only  dependent 
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on  gas  temperature  but  is  also  a  function  of  gas  pressure.  Estimations  based  on  the 
assumption  that  the  cooling  water  is  not  boiling  in  the  channel  give  a  ranee  of  regulated 
helium  thermal  conductivities  between  0.073-0.087  W/mK  at  a  gap  pressure  of  10  torr. 
This  range  is  approximately  30%  of  the  thermal  conductivity  value  calculated  assuming  a 
temperature  dependence  only.  The  regulated  helium  gap  is  one  of  the  most  critical 
resistances  in  the  thermal  network. 

The  primary  mode  of  heat  transfer  into  the  water  channel  is  forced  convection  heat 
transfer.  Determining  the  amount  of  heat  transferred  into  the  water  by  forced  convection 
requires  calculating  the  convective  heat  transfer  coefficient.  This  coefficient  is  dependent 
on  the  flow  characteristics  in  the  channel,  the  cooling  channel  geometry  and  the 
thermophysical  properties  of  the  cooling  water.  The  Reynolds  number  is  calculated  to  be 
1480  which  indicates  that  the  flow  rate  is  in  the  combined  laminar  flow  regime.  This  is  a 
very  difficult  regime  in  which  to  accurately  calculate  the  heat  transfer  coefficient  due  to 
the  importance  of  the  buoyant,  conductive  and  convective  forces  in  the  channel.  Based  on 
existing  correlations,  the  convective  heat  transfer  coefficient  is  determined  to  fall  between 
the  range  of  528-1012  W/m2K. 

Another  critical  resistance  in  the  thermal  network  is  the  cesium  interelectrode  gap. 
Conduction  heat  transfer  occurs  through  the  cesium  vapor  and  can  be  defined  in  terms  of 
the  cesium  vapor  pressure,  the  gap  size  and  the  surface  temperatures  of  the  emitter  and 
collector.  Due  to  the  high  emitter  temperatures,  radiation  heat  transfer  is  significant 
through  this  gap.  The  gap  heat  transfer  characteristics  are  directly  affected  by  the 
thermionic  process.  Coupling  the  thermal  model  with  a  thermionic  model  is  crucial  to 
estimating  emitter  temperatures  based  on  data  collected  from  the  test  stand. 

The  TIT  AM  computer  code  provides  a  means  of  simulating  steady-state  and 
transient  operation  of  a  TOPAZ-II  TFE  in  the  reactor  core.    TIT  AM  can  be  utilized  to 
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estimate  the  operating  characteristics  of  a  TFE  in  the  test  stand  by  inputting  the  active 
zone  power,  average  collector  sleeve  temperature  (corresponding  to  the  NaK 
temperature),  and  an  assumed  NaK  flow  rate  (based  on  typical  values).  The  collector 
sleeve  temperature  calculated  by  TITAM  is  766  K  compared  to  a  temperature  of  764  K 
calculated  using  the  thermal  network.  These  temperatures  calculations  are  based  on 
simple  heat  transfer  principles  which  provides  good  correlation  between  the  two 
calculations.  TITAM  can  also  be  used  to  estimate  the  emitter  surface  temperature  based 
on  the  input  data,  however,  the  thermal  resistance  of  the  cesium  interelectrode  gap  must 
be  defined  before  the  thermal  analysis  presented  in  this  thesis  can  be  used  to  provide  an 
estimate  of  the  emitter  surface  temperature  to  compare  it  to. 

Differences  in  the  TITAM  output  and  the  thermal  analysis  results  can  be  attributed 
to  several  reasons.  The  most  important  of  these  is  that  TITAM  simulates  a  TFE  operating 
in  the  TOPAZ-II  reactor  core  and  not  in  the  test  stand.  Several  other  differences  include 
cesium  reservoir  temperature,  TFE  load  voltage,  TFE  load  current,  and  TFE  efficiency. 
To  understand  the  discrepancies  between  the  two  results  requires  a  detailed  examination 
of  the  TITAM  code  and  integration  of  the  thermal  model  developed  in  this  thesis  with  the 
thermionic  model  of  the  code. 
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VII.  RECOMMENDATIONS 

The  TOPAZ-II  TFE  test  stand  provides  an  excellent  vehicle  for  investigating  the 
thermionic  characteristics  of  a  single  TFE;  however,  understanding  the  thermal 
characteristics  of  the  test  stand  and  how  its  thermal  behavior  affects  the  thermionic 
process  is  crucial  in  being  able  to  correlate  test  stand  results  to  the  entire  reactor  core. 
Several  improvements  to  the  test  stand  and  additional  testing  could  greatly  aid  in 
improving  the  thermal  model  developed  in  this  thesis. 

A  detail  investigation  of  the  flow  characteristics  of  the  cooling  water  would 
provide  a  better  estimation  of  the  convective  heat  transfer  coefficient.  This  could  include 
moving  the  inlet  and  outlet  cooling  water  thermocouples  inside  the  vacuum  chamber 
closer  to  the  test  sections  water  jacket.  Using  more  accurate  temperature  measurements 
would  allow  increasing  the  water  flow  rate  into  a  better  understood  and  documented  flow 
regime.  Testing  of  the  cooling  water  characteristics  could  be  accomplished  using  a  known 
power  input  into  the  degassing  heaters  and  would  not  require  the  operation  of  the  TFE. 
Thermocouple  instrumentation  inside  the  water  channel  itself,  would  also  prove  useful. 

Further  understanding  the  pressure  dependence  of  the  regulated  helium  gap's 
thermal  conductivity  is  critical  to  the  accuracy  of  the  thermal  model.  Once  the  thermal 
characteristics  of  the  regulated  helium  gap  and  the  cooling  water  channel  are  well  defined, 
a  more  accurate  heat  balance  can  be  performed  between  the  collector  sleeve  and  the 
cooling  water. 

Exploring  the  end  effects  at  the  upper  and  lower  helium  chambers  would  provide  a 
better  understanding  of  the  2-dimensionality  of  the  thermal  system.  Experiments  could 
easily  be  performed  without  the  operation  of  the  TFE  to  examine  how  these  end  chambers 
are  affecting  the  axial  temperature  distribution  along  the  collector  sleeve.    Operating  the 
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heaters  for  one  or  both  of  these  chambers  during  operation  of  the  TFE  provides  an 
additional  heat  input  into  the  thermal  model;  therefore,  the  input  power  to  these  heaters 
should  be  recorded. 

Experimental  information  is  only  valuable  if  the  accuracies  of  the  instrumentation 
and  data  acquisition  system  are  completely  defined.  The  integration  of  Russian  and 
American  equipment  is  inherent  to  the  TOPAZ  project  and  increases  the  difficulty  in 
developing  an  accurate  mode!  for  the  test  stand.  Improved  documentation  of  the  test 
stand  measurement  errors  should  accompany  any  future  experimental  work. 

Experimental  data  should  be  collected  from  the  test  stand  at  a  predetermined  TFE 
output  voltage.  This  voltage  level  should  be  based  on  supplying  a  specified  bus  voltage 
(i.e.  28  VDC  or  30  VDC).  This  must  be  a  consideration  since  operating  the  TFE  at  the 
maximum  output  power  (or  maximum  efficiency)  might  not  produce  useful  voltage  to 
supply  a  power  distribution  bus. 

TITAM  is  a  valuable  tool  that  should  be  modified  further  for  use  with  the  test 
stand.  Utilizing  the  thermionic  emission  model  from  TITAM  along  with  the  thermal 
model  developed  in  this  thesis  can  provide  a  detailed  understanding  of  the  thermal 
characteristics  of  the  test  stand.  This  would  allow  for  a  more  direct  correlation  of  the  test 
stand  results  with  the  operation  of  the  overall  TOPAZ-II  reactor  system. 
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APPENDIX  A.  THERMAL  PROPERTIES 


The  thermal  conductivities  and  emissivities  of  various  test  stand  components  are 
graphed  and  are  defined  in  terms  of  temperature.  Most  of  the  information  was  available  in 
Paramonov  and  El-Genk,  1994,  but  some  was  derived  from  reference  book  data  tables. 
References  are  provided  in  the  titles  of  each  figure.  A  polynomial  curve  fit  of  the 
tabulated  data  accomplished  using  Mathcad Plus  5.0  (MathSoft,  Inc.,  1994). 
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APPENDIX  B.  EXPERIMENTAL  DATA  TABLES 


Data  were  collected  in  support  of  this  thesis  during  two  separate  experimental 
sessions.  The  first  occurred  between  15-26  March,  1994,  when  data  were  collected  over  a 
TISA  heater  input  power  range  of  1000  W  to  3600  W.  The  primary  purpose  of  the  runs 
completed  during  this  period  was  to  establish  the  operating  characteristics  of  the  test  stand 
and  its  associated  TFE  from  normal  operating  range  to  very  low  input  power  levels.  This 
information  is  being  used  in  a  follow-on  thesis  to  examine  the  thermionic  characteristic  of 
a  TFE  in  the  test  stand. 

The  second  experimental  session  occurred  between  14-18  August,  1994  Prior  to 
this  session,  the  entire  test  stand  data  acquisition  system  was  converted  from  a  computer 
hardware  and  software  system  that  was  furnished  by  the  Russians  to  a  more  contemporary 
system  that  was  designed  and  assembled  by  the  US,  TIP  personnel.  This  system  update 
provides  for  much  easier  data  acquisition  and  test  stand  control.  Although  the  data 
acquisition  system  changed  between  the  two  experimental  procedures,  the  test  stand 
operating  conditions  were  maintained  as  close  as  possible  to  allow  for  data  correlation 
between  the  two  collection  periods.  The  primary  purpose  of  the  second  experimental 
procedure  was  to  support  thermal  model  calculations  by  investigating  the  relationship 
between  regulated  helium  gap  pressure  and  collector  sleeve  temperature  and  to  collect 
data  a  various  power  levels  for  an  optimized  cooling  water  flow  rate.  Relevant  data  to  the 
thermal  analysis  are  tabulated  in  this  appendix. 

The  appendix  is  separated  into  three  time  periods,  each  time  period  consisting  of 
three  tables.  The  shaded  columns  contain  data  that  were  recorded  from  the  test  stand  and 
the  unshaded  columns  contain  derived  data.  All  of  the  data  are  collected  at  a  cesium 
reservoir  temperature  of  approximately  255  °C  (around  0.5  torr  of  pressure). 
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The  second  column  of  the  first  table  is  TIS  A  active  zone  power  and  is  determined 
by  multiplying  the  TISA  input  power  by  88%.  The  fifth  column  is  TFE  output  power  and 
is  determined  by  multiplying  the  TFE  output  voltage  by  its  output  current.  The  column 
indicated  as  column  "A"  is  the  power  loss  from  the  TFE.  These  values  are  found  by 
subtracting  the  TFE  output  power  from  the  TISA  active  zone  power.  TFE  efficiency  in 
column  seven  is  the  result  of  dividing  the  TFE  output  power  by  the  TISA  active  zone 
power  and  multiplying  by  100. 

Column  eight  of  the  first  table  is  the  average  collector  sleeve  temperature  and  is 
the  same  value  as  in  column  five  of  the  third  table.  Performing  a  heat  balance  on  the 
cooling  water  provides  the  data  in  the  columns  indicated  as  column  "B"  and  "C".  There  is 
one  column  for  each  of  the  two  output  thermocouples.  These  values  are  the  same  values 
as  given  in  columns  12  and  13  of  the  third  table.  The  %  error  in  columns  1 1  and  13  of  the 
first  table  are  the  percent  difference  between  the  cooling  water  heat  balance  (columns  "B" 
and  "C")  and  the  power  balance  (column  "A") 

The  second  table  presents  the  11  collector  sleeve  thermocouple  measurements. 
See  Figure  5-1  for  the  relative  axial  location  of  the  thermocouples  along  the  collector 
sleeve.  Columns  2  through  6  of  the  third  table  provide  relevant  information  concerning 
the  thermocouple  measurements.  Column  five  gives  the  average  of  all  1 1  collector  sleeve 
thermocouples  and  column  six  gives  the  average  of  thermocouples  T57  through  T66.  This 
excludes  the  end  point  temperature  measurements  of  T56  and  T67. 

Column  eight  of  the  third  table  contains  the  measured  water  flow  rate  (in  gpm) 
converted  to  a  cooling  water  mass  flow  rate  (in  kg/s).  Columns  12  and  13  of  the  third 
table  indicate  the  results  of  the  cooling  water  heat  balance  (utilizing  Equation  5.10). 
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The  tables  dated  13-14  August  1994  are  taken  over  a  TISA  input  power  of 
1000-3600  W  At  each  power  level  the  regulated  helium  gap  pressure  was  adjusted  to  its 
optimum  value.  This  maintains  a  relatively  constant  collector  sleeve  temperature  over  the 
entire  power  range.  The  high  %  error  for  the  1000  W  data  run  is  due  to  the  small 
temperature  difference  between  the  cooling  water  inlet  and  outlet.  This  indicates  how 
sensitive  the  cooling  water  heat  balance  is  to  this  temperature  difference. 

The  tables  dated  14-16  August  1994  provides  data  collected  over  three  power 
ranges  with  varying  regulated  helium  gap  pressure  at  each  power  level.  The  first  table 
shows  that  for  a  constant  input  power  level,  the  collector  sleeve  temperature  decreases  as 
the  regulated  helium  gap  pressure  increases.  This  proves  how  important  understanding 
the  pressure  dependence  of  the  regulated  helium  gap  thermal  conductivity  is  to  the  thermal 
analysis  of  the  test  stand. 

The  tables  dated  17-18  August  1994  provides  a  set  of  data  runs  under  similar 
conditions  as  presented  in  the  first  set  of  tables.  A  complete  set  of  all  data  acquired  during 
both  experimental  periods  is  archived  at  the  TIP  facility  in  Albuquerque,  New  Mexico. 
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APPENDIX  C.  THERMAL  NETWORK  CALCULATIONS 


The  following  calculations  and  the  calculations  in  Appendix  D  were  completed  using 
Mathcad  Plus  5.0  (Mathsoft,  Inc.  1994).  This  calculational  software  uses  a  colon  followed 
by  an  equal  sign  (:=)  to  mean  "defined  by"  vice  an  equal  sign  (=)  to  mean  "equal  to"    All 
variables  must  be  defined  prior  to  being  used  in  an  equation.  After  a  variable  is  defined,  an 
equal  sign  can  be  used  to  determine  the  variables  value. 

A.      THERMAL  MODEL  INPUTS 

1.  INPUT  POWER  TO  THE  TISA  HEATER 

QTISA  r  3412  watt 

2.  POWER  TO  THE  ACTIVE  ZONE 

Qaz  =  QtISA'0-88         (2^=3.003. 103    -watt 

3.  TFE  OUTPUT  POWER 

Pout  =  108.9  watt 

4.  POWER  LOSS  FROM  TFE 

Q  loss  -  Q  az  -  p  out       Q  loss  =  2  894*  1q3    >watt 

5.  RADIAL  DIMENSIONS  FOR  THE  TEST  STAND 

r\   --  0.0098m         r4  =  0.0119m       r7  =  0.01495m       r  10  =  0.0300m 

r2  =  0.0103  m         r5  =  0.01195  m      rg  =  0.0150  m         rn   =  0.0325  m 
r3   =  0.01175  m       r6  =  0.01370  m     r9  =  0.0295  m         r  12  =  0.0375  m 

6.  LENGTH  OF  THE  TEST  SECTION 

z  :=  0.3  m 

7.  DIMENSIONAL  CONSTANT  FOR  CALCULATIONS 

c  :=  2-tc-z        c  =  1.885-m 
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8.  COLLECTOR  SLEEVE  AXIAL  TEMPERATURE  DISTRIBUTION 

T56   =  424.49  K         T  59   =  483.9  K  T  62   =  433.73  K       T  66  =  362. 95K 

T57  =  515.64  K         T60  =  482.57  K         T  63    =  439.32  K       T  67  =  354.29  K 
T58   =  498.46  K         T61   =  490.12  K         T  65   =  376.42  K 

9.  AVERAGE  COLLECTOR  SLEEVE  TEMPERATURE 

The  average  collector  sleeve  temperature  is  calculated  without  using  the  end  point 
temperatures,  T(56)  and  T(67). 


T6c 


57  +  l  58  +  l  59  +  l  60  +  l  61  +  l  62  +  l  63  +  l  65  +  l  66 


9 
T6  :=  T6c  ♦■  273.15  K  T6=726.829-K 

B.      FROM  THE  COLLECTOR  SLEEVE  OUTWARD 

1.  STAINLESS  STEEL  COLLECTOR  SLEEVE 

a.  Thermal  Conductivity  of  Stainless  Steel 

i  ,~«/wat*\     ~~,  /  watt  \  _  ,  .««,«  watt 

kss  =  ,2.0-(--).0.0l|_).T6         kss  =  ,9.268.-- 

b.  Temperature  Drop  Across  Collector  Sleeve 

M 

Qlossln  — 

\r6/ 

AT  ss  = —  AT  ss  =  6  957*K 


ck 


SS 


c.  Temperature  at  node  7 


T6-ATSS  T7=719.872-K 
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2.  UNREGULATED  HELIUM  GAP 


a.  Thermal  Conductivity  of  a  Helium  Continuum 

i                  o  in  a  ia  4  /  watt  \  x  o  ii  in"7  /  watt  \  t     2 

kHel    =9.376-10        -T7-  9.1310    ■-— -)-T7 

mlC/  \m-K 


"h*       ;.5.10->".M.T73-, .264.10"". pi- T/ 
m-K  /  \mK 


«  «,  ,«-2  /watt\      ,  ,  ,  ^  „ watt 

kHe      -9.0610  M--NkHcrkHe2       kHe  =0.282  —  - 
\m  K/  m  K 


b.  Temperature  Drop  Across  the  Unregulated  Helium  Gap 

/r8 
Qlossln  - 

ATHp  = i-^- 


c.  Temperature  at  Node  8 


T8  =  T7  -  ATHe  T8  =701.725-K 


3.  COPPER  CONDUCTOR 


a.  Thermal  Conductivity  of  Copper 


,,%„  „„^  watt      ^  -,, ,~~  ,„-2  /  watt  \  ^         ,  „„„  ,„  watt 

kCu  =422.286  —--6.71429  10 -T  8      kCu  =  375.17  —  - 

mK  mK2  mK 
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b.  Temperature  Drop  Across  the  Copper  Conductor 

/r9 
Qlossln|  — 

ATCu       -r^  ATCu  =  27 

c  kCu 

c.  Temperature  at  Node  9 

T9=T8-ATCu  T9=  698.958-K 

4.  REGULATED  HELIUM  GAP 


The  thermal  conductivity  of  the  helium  in  the  regulated  gap  is  first  estimated  using 
the  equation  for  a  continuum  gas. 

a.  Thermal  Conductivity  of  a  Helium  Continuum 


4  /  watt  \  _        «  ,-  ,«-7  /  watt  \  _    2 
' T9  -  9.13  10       ;  ' 

mK2  \rn-K' 


He3   =9.37610-        —    T9-  9.13  10'-  -]-T9 


kHe4  =  5.510  10p^    T93  -  1.26410  ,3    ^  - 

\m-K4/  \m-K5 

„  ~s  ,~-i  /watt\     ,  ,  ,  «  „„,  watt 

kHereg  '  "906  10    \~\  *  kHe3  ♦  kHe4       kHereg  =0  216'~^ 


b.  Temperature  Drop  Across  the  Regulated  Helium  Gap 

'rio\ 


Qloss  ln    | 

AT  Hereg  =  —  j— ^  ^  Hereg  =  93.366-K 

CK  Hereg 
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c.  Temperature  at  Node  10  Based  on  a  Helium  Continuum  in  the  Gap 

T  lOcont      T  9  -  AT  Hereg     T  10cont  =  605. 592 -K 

This  temperature  will  lead  to  a  water  channel  inner  wall  temperature  of  around 
325  °C.   Since  this  is  approximately  200  °C  higher  than  the  boiling  point  of  water,  the 
helium  thermal  conductivity  must  be  much  lower  than  its  continuum  value.  From  Figure 
5-5,  the  helium  thermal  conductivity  is  estimated  to  be  0  08  W/mK  for  a  wall  temperature 
of  80  °C.  Using  a  temperature  range  of  80  +  30  °C  gives  a  low  thermal  conductivity  of 
0.073  W/mK  and  a  high  thermal  conductivity  of  0.087  W/mK.  The  temperature  drop 
across  the  regulated  helium  gap  is  recalculated  based  on  the  low,  estimated,  and  high 
thermal  conductivities. 

d.  Estimated  Temperature  Drop  Across  the  Regulated  Helium  Gap 


10 


n  A„„  watt  ,oss 

"Helow  ■  0.073--        &r ^  .         ^2 

KHelow 
rl 


Qlossln 


Heest 


0.08 


watt 
m-K 


AT 


Heest 


ck 

Qlossln 


Heest 
r10 


k  Hehigh  -0.087- 


m 


K     "  Hehi8h      —  k 


Hehigh 


ATHelow=353441-K 


ATHeest=322515-K 


AT  Hehigh  =296.565.K 


e.  Estimated  Temperature  at  Node  10 


T  lOlow  =  T9  "  ATHelow  T  lOlow  "345.517-K 


T10  -  T9-  ATHeest 


T10  -376.443-K 


T  lOhigh  -  T  9     AT  Hehigh  T  lOhigh  =  402.393 -K 
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5.  Stainless  Steel  Water  Jacket 


a.  Thermal  Conductivity  of  Stainless  Steel 


k^lM«-O01 


watt 


mK' 


T 


lOlow 


wjlow 


=  15.455 


watt 
mK 


wj 


,~  ~  /watt\      n  n,  I  watt  . 
--  12.0 +  0.01  I I 


mK, 


mK' 


10 


*r    -,sa      Watt 

kwj  -15-764.-- 


kwjhigh  =  120, 


watt  \     „  n ,  /  watt  \  _  ,  , ,  «„ ,  watt 

+  0.0l( ^)T10high        kwjhigh  =  16.024 

mK 


mK 


b.  Temperature  Drop  Across  the  Water  Jacket 

rll 


QiossH 


AT 

wjlow 


10, 


ck 


Qlossln 


wjlow 

rll 


AT 


10; 


wj 


ck 


wj 


AT  wjlow  =795"K 


AT    •  =7.795-K 


Qlossln 


11 


AT 


wjhigh 


T\ol 


ck 


wjhigh 


AT  wjhigh  =  7668-K 


c.  Temperature  at  Node  11 


T  1  How  =  T  lOlow  "  AT  wjlow  T  1  How  "  337.567-k 


11  ;=  l  10 


ATwj 


Tn  =368.649-K 


T  llhigh  ^  T  lOhigh  "  AT  wjhigh        T  1  lhigh  =  394.724-K 
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These  calculations  yield  a  water  channel  inner  wall  temperature  of  between 
64-122  °C.   Using  the  Dittus-Boelter  correlation  convective  heat  transfer  coefficient  of 
1012  W/m2K  as  the  upper  limit  and  the  Seider-Tate  correlation  convective  heat  transfer 
coefficient  of  528  W/nrK  as  the  lower  limit,  the  last  resistance  in  the  network  can  be 
calculated  and  a  bulk  cooling  water  temperature  can  be  calculated. 

6.  WATER  CHANNEL 

a.  Convective  Heat  Transfer  Coefficient  and  Heat  Transfer  Area 
h,„„,  .  528-^1  A  ht,ow  .0.07  m2 


low 


m2 "K  Ahthigh  =  0.125  m2 


watt 

hhigh  =  1012^—  A  htlow  *  A  hthigh  rt  rtrt£1     - 

*  m2K  Aht  = 5-       Aht=0.098.m 

b.  Temperature  Drop  from  the  Wall  to  the  Bulk  Fluid  Temperature 

Qloss 

ATwlow^T-r—  ATwlow  =  56209'K 

Aht  hlow 


Qloss 
hthhigh 


ATwhigh  =  ^-.^77  ATwhigh  = 


c.  Bulk  Water  Temperature 

Tbulklowl    =  Tlllow  "  ATwlow  Tbulklowl  =281.357-K 

Tbulklow2  =  T  lllow  '  ATwhigh  Tbulklow2  =308.24-K 

Tbulkl   =  Tll  "  ATwlow  Tbulkl  =312.439-K 

T  bulk2  =  T  1 1  "  AT  whigh  T  bulk2  =  339.322-K 

Tbulkhighl   =  T  11  high  _  ATwlow  Tbulkhighl  =  338.5 15-K 

Tbulkhigh2  =  T  11  high  "  AT  whigh  Tbulkhigh2  =  365  39S'K 

These  calculations  give  a  bulk  cooling  water  temperature  range  of  8-92  °C. 
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C.      FROM  THE  COLLECTOR  SLEEVE  INWARD 

1.  STAINLESS  STEEL  COLLECTOR  SLEEVE 
a.  Thermal  Conductivity  or  Stainless  Steel 


,~  ~  /watt\      «  „,  /  watt  \  „  ,  .««,,„  watt 

kss.,2.0.(--).0.01.—    T6         lcss -19268.-- 


b.  Temperature  Drop  Across  Collector  Sleeve 


,r6 

Qlossln|  — 
r5 


AT  ssi  = r^~^         AT  ssi =  10  888<K 

ckSS 
c.  Temperature  at  node  5 


T  5  ■=  T  6  +  AT  ssi  T  5  =  737.717-K 


2.  UNREGULATED  HELIUM  GAP 


a.  Thermal  Conductivity  of  a  Helium  Continuum 


n  -.-.^  ,«-4  /  watt  \  _        .  ,„      -7  /  watt  \  _    2 
Hel.  =  9.376  10 T5-  9.13  10       --T52 

\nvK  /  \m-K  / 


kHe2i  ■.  s.3-10-10-p£   T53  -  1.264- 10- "p^V 
\m-K  /  \m-K 


Hei  =  "906  10  2 •[—)  +  ^Heli  +  kHe2i    kHei  =0.288-^ 

\mK/  m  K 
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b.  Temperature  Drop  Across  the  Unregulated  Helium  Gap 

Qioss1"!  — 

*THei         -ZV^         ATHei=2238: 
c  KHei 

c.  Temperature  at  Node  8 

T  4   =  T  5  +  AT  Hej  T  4  =  760.099-K 

3.  ALUMINA  INSULATION 


a.  Thermal  Conductivity  of  Stainless  Steel 

-,  r,  W^tt     rt  ,__   watt  _        ,  „  ,A-4   watt  „,    2 
kaluml   =74.9--  -0.178 -T4  +  1.7510    -    — T4 

mK  mK2  mK3 

kalum2=-7.7810-8^.T4'.,.3,..0"^L.T; 

m-K  m-K 

,  watt 

kalum  =  kaluml  +  k  alum2  kalum~10916 —  ~ 

m  K 


b.  Temperature  Drop  Across  the  Alumina  Insulation 

h 

Qloss1"!  — 
r3 


AT  alum  ■  —  7-^^  AT  alum  =  '  784'K 

c  Kalum 


c.  Temperature  at  node  3 


T3   =T4^Talum  T3=761.883.K 
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4.  COLLECTOR 

a.  Thermal  Conductivity  of  Polycrystalline  Molybdenum 


,,«  watt      ,  _  ,„-2   watt   ^  ,  „„  „„,,  watt 

kpolyMo      >W~j      15  1° J-T3         kpolyMo  = 98  572~ ^ 

m*Jv 

b.  Temperature  Drop  Across  the  Collector 

Qloss1"  — 

ATcoU--7  ATCO||=2.051.K 

c  K  polyMo 


c.  Temperature  at  Node  2 


T2  =  T3  +  ATcoll  T2=763.934.K 
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APPENDIX  D.  WATER  CHANNEL  CALCULATIONS 

WATER  CHANNEL  HEAT  TRANSFER  AREA 

1.  Dimensional  Constants 

rji    =  0.0325  m  Inner  Cylinder  Radius  in  m 

r  13   =  0.0425  m  Outer  Cylinder  Radius  in  m 

z  =  0.326  m  Water  Jacket  Length  in  m 

rc  =  0.005  m  Radius  of  the  Flow  Guidance  Coil 

d  cj  =  0.075m  Test  Stand  Diameter  to  center  of  the  Water 

Channel  in  m 
9  =  13  deg  Cant  Angle  of  Flow  Guidance  Coil  in  degrees 

N  =  6.5  Number  of  Revolutions  of  the  Flow  Guidance 

Coil 

2.  Inner  Cylinder  Surface  Area  in  m2: 
A j  =  2  7t  r  ji  z 

Aj  =  0.067-m2 

3.  Outer  Cylinder  Surface  Area  in  m2: 
A0  =  2-7C-r13  z 

A0  =  0.087-m2 

4.  Surface  Area  of  the  Flow  Guidance  Coil  in  m2: 

dcl 
d  eff  -  Effective  Test  Stand  Diameter  based  on  the  Flow 

cos(9)  Guidance  Coil  Cant  Angle  in  m 

deff  =  0.077- m 

L  :=  N  n  d  eg»  Length  of  the  Flow  Guidance  Coil  in  m 

L  =  1.572-m 

Cc  =  2  n  t  c  Circumference  of  the  Flow  Guidance  Coil  in  m 

Cc  =0.031-m 

Ac  ;=  L  Cc  Surface  Area  of  the  Flow  Guidance  Coil  in  m2 

Ac  =  0.049- m2 
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5.  Total  Surface  Area  in  Water  Channel  in  m2: 


At  .-  Aji- A0*  Ac 

At  =0  203-m2 


B.      HEAT  TRANSFER  CORRELATIONS 


1.  Definitions  of  Physical  Properties  of  Water  at  15°C  (Welty  et  al.,  1984): 

p  =  1000.0    —  Density 

\m3/ 

c  „  =  4185.0 Specific  Heat  Capacity 

P  \kg-K/ 

1 1  6  10  4  I     —     1  Dynamic  Viscosity 

Im-sec/ 

/    2\  . 

116  10"7 ■  —  Kinematic  Viscosity 

\sec/ 

k  =  0.59 Thermal  Conductivity 

\m-K/ 

Pr  ■=  — ^—  Pr  =  8.228       Prandtl  Number 

k 

2.  Total  Cooling  Water  Flow  Area  in  m2: 

Af  =  (0.035  m)  (0.010  m)  -  k  (0.005  m)2 
Af  =  2.715-10"4    -m2 

3.  Wetted  Perimeter  in  m: 

Pw  =  (2  0.035  m)  +  (2-tct 


w  V~    v,.^^^    „V  ,  s.    „   .    c 


Pw=0.101-m 
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4.  Equivalent  Hydraulic  Diameter  in  m: 


De  =4 


Pw 


De  =  0.011-m 


5.  Cooling  Water  Flow  Velocity  in  m/s: 

-^-1        V  =  4.353- 10  5    -m3 

Mnl  Flow  rate 


V  =  0.69-1  — I        V  =  4.353-10  5    -m3 -sec"1         - —  Measured  Volumetric 


V 

v  =  — 

Af 

v  =  0.16-m-sec 
6.  Reynolds  Number: 
Devp 


Re 


3 


Re  =  1.48-10 

7.  Seider-Tate  Relation  for  Laminar  Flow: 

Pe  ■=  Re-Pr         Pe  =  1.218-104       Peclet  Number 

\i       --  352  10  6  (  —  ^-  |  Wall  Temperature  Fluid  Viscosity  for  80°C : 

\m-sec/ 

\_ 

k  I      De\3'  ■■   \°-14 


hST 

hST  =528.245- 


De""\"    L/     [»„) 
watt 


m2K 
8.   Dittus-Boelter  Correlation  for  Fully  Developed  Turbulent  Flow: 

.  k  0.8  D  0.4 

hnn  = 0.023  Re      Pr 

D  a 
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